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Abstract. We set up an algebraic theory of multivariable integration, based on a hierar¬ 
chy of Rota-Baxter operators and an action of the matrix monoid as linear substitutions. 

Given a suitable coefficient domain with a bialgebra structure, this allows us to build an 
operator ring that acts naturally on the given Rota-Baxter hierarchy. We conjecture that 
the operator relations are a noncommutative Grdbner basis for the ideal they generate. 

1. Introduction 

1.1. Motivation. The notion of integral operator plays a fundamental role in Analysis, 
Physics and Stochastics. Apart from some early beginnings, the systematic study of in¬ 
tegral operators and integral equations started in earnest by the end of the nineteenth 
century, at the hands of towering hgures like Hilbert, Fredholm and Volterra. 

Despite their indubitable origin in Analysis, it is often prohtable to study integral op- 
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algebras [17]. The situation is similar to the case of differential algebra, dehned to be an 
(associative) algebra equipped with a derivation. Its study originated from the algebraic 
study of differential equations by Ritt [31, 32] starting in the 1930s and has since developed 
into a vast area of mathematics research with broad applications, including the mechanical 
proof of geometric theorems [25, 44]. An algebraic structure that encodes integration hrst 
appeared as a special case of a Rota-Baxter algebra, dehned to be an algebra R equipped 
with a linear operator P on R satisfying the Rota-Baxter axiom 

P{u)P{v) = P{uP{v)) -I- P{P{u)v) + XP{uv) 

for all u,v ^ R. Here A G R is a prehxed constant called the weight, which is zero in the 
case P = but nonzero for its discrete analog (partial summation). Let us now name a 
few important areas where this approach proved to be successful: 

• The original formulation of Rota-Baxter algebra (hrst called Baxter algebras) was 
in the context of fluctuation theory, when Glen Baxter gave a new proof [4] of 
the so-called Spitzer identity, based on the Rota-Baxter axiom. This line of re¬ 
search continues to this day; see the references [12, 61, 121, 163, 174] in the survey 
monograph [17]. 

• The use of Rota-Baxter algebra for problems in combinatorics has been hrmly 
established and widely popularized [38, 39] by Gian-Carlo Rota (which is why they 
are partly named after him). For a typical example, see [17, §3.3] and the references 
therein. 

• One of the main applications of Rota-Baxter algebras is the operator form of the so- 
called classical Yang-Baxter eguation, named after the physicists Ghen-Ning Yang 
and Rodney Baxter (not Glen Baxter). For an overview and current developments 
in this area, see the references [14, 26, 46, 127, 164] in [17]. 

• Another important application in Physics is in the renormalization of perturbative 
quantum held theory, in particular the algebraic Birkhoff decomposition established 
by Connes and Kreimer [11]. A survey of current literature in this domain is given 
in [17, §2.5]. 

• The relation of Rota-Baxter algebras and multiple zeta values and their general¬ 
izations is a fascinating topic [20, 21, 45]. For more on this topic, the reader may 
read [17, §3.2] and the literature survey in [17, §3.4]. 

• Algorithmic tools for solving and factoring linear boundary problems [13, 40] are 
developed [33, 35, 37] based on integro-differential algebras (Rota-Baxter algebras 
with a suitable differential structure). The solution of the boundary problem is 
denoted by the Green’s operator, an element in a ring of integral operators. 

In this paper, we shall employ Rota-Baxter algebras for creating an algebraic theory of 
multivariable integation under the action of linear substitutions. As far as we are aware, 
this is the hrst time that integral operators and substitutions are studied from a systematic 
algebraic perspective. 
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Our original motivation comes from the last of the application areas listed above: One 
may build up an abstract theory of linear boundary problems and their Green’s opera¬ 
tors [30], “abstract” in the sense that it abstracts from the ring-theoretic structures (deriva¬ 
tions, Rota-Baxter operators, characters), retaining only their linear skeleton. The advan¬ 
tage of the abstract theory is that it encompasses both linear ordinary differential equations 
(LODEs) and linear partial differential equations (LPDEs), both scalar and vector prob¬ 
lems, both continuos (differential) and discrete (difference) equations. However, unlike the 
integro-differential setting of [35], it does not—and cannot—provide algorithmic tools. 

For developing an algebraic theory of boundary problems for LPDEs, along with the 
required algorithmic tools, it is hrst of all necessary to develop a suitable ring of par¬ 
tial integro-differential operators capable of expressing the Green’s operators of (some of) 
these boundary problems. An adhoc example of such an operator ring—together with 
applications to solving certain simple LPDE boundary problems—can be found in [34]. 

It must be emphasized that a suitable ring of partial integro-differential operators is not 
just the tensor product of the corresponding ordinary integro-differential operator rings. 
Such a ring would be too weak for expressing the Green’s operators of even very simple 
LPDEs like the inhomogeneous wave equation uu — Uxx = / mentioned in [30, §7], and the 
problems treated in [34]. The reason for this is evident from a geometric perspective: The 
tensor product provides only axis-parallel integrators while most Green’s operators require 
integration over oblique domains (slanted lines, triangular areas, skewed volumes). This 
contrasts with differential operators, where any directional derivative can be represented— 
via the chain rule—as a linear combination of partial (“axis-parallel”) derivatives. 

For rectifying this defect, it is sufficient to add linear substitutions to the operator ring. 
Their algebraic interation with the axis-parallel integrators is regulated by (a special case 
of) the substitution rule of integration. The precise algebraic formulation turns out to 
be surprisingly delicate, and the current axiomatization (Dehnition 2.3) should be seen 
as a hrst attempt in this direction. Our central question can be stated thus: “What is 
the algebraic structure of linear substitutions under partial integral operators?” See the 
Gonclusion for some further ideas on how to improve and generalize the current setup. 

The presence of linear substitutions also makes the choice of the coefficient algebra much 
more delicate (see the beginning of Section 3 for a short explanation). To incorporate 
closure under linear substitutions into the algebraic setting, we require the coefficients to 
form a Hopf algebra with a suitably compatible Rota-Baxter operator; we call the resulting 
structure a Rota-Baxter bialgebra (Dehnition 3.5). 

For simplifying the further treatment, we have decided to keep diherential operators 
out of the picture. In other words, we are building up rings of partial integral operators 
(with linear substitutions and coefficient functions), rather than partial integro-differential 
operators. In view of the applications for boundary problems, this is clearly not sufficient 
but it provides a crucial intermediate stage. In fact, the adjunction of partial diherential 
operators (coming from derivations of the base algebra) leads to new operator relations 
that are much simpler than the ones treated here: Loosely speaking, diherential operators 
“march through” to the far right (unlike integral operators) since the chain rule leads 
to a full decoupling, as indicated by our earlier remarks above. We shall develop the 


4 


MARKUS ROSENKRANZ, XING GAO, AND LI GUO 


resulting operator ring in a subsequent paper, along with its application to expressing 
Green’s operators of LPDE boundary problems. 

Apart from this deliberate restriction, the operator rings constructed in this paper are 
a vast generalization of the adhoc ring created in [34, 36] for the purpose of illustrating 
certain techniques for solving LPDE boundary problems. To wit, the current ring is not 
based on the specihc model of multivariate analytic functions as in [34], and it allows 
for a wide class of coefficient algebras (arbitrary Rota-Baxter bialgebras) rather than just 
constant coefficients. 

From the viewpoint of Analysis, one may also pose the following question: Let R be the 
subring of the ring of linear operators on the complex function space C'°°(]R"') generated by 
the partial integral operators, the multiplication operators induced by certain coefficient 
functions, and all linear substitutions. If P denotes the set of all these generators, there 
is an evident map C(P) —)■ R that “evaluates” the indeterminates by the corresponding 
generators. The kernel TZ of this map encodes the ideal of operator relations satished in 
Analysis. If these are known, one has a full algebraization in the sense that R = C(P)/7^. 
While we have not reached this goal in the present paper, the relations TZ' used in our 
formulation of the operator ring (Dehnition 4.8) are clearly an important subset of TZ, and 
there are good reasons to conjecture that actually TZ = TZ'. However, this will presumably 
need an analytic proof that has to be approached elsewhere. 


1.2. Structure of the Paper. In Section 2 we introduce the main structures needed 
for the algebraic description of multivariable integration and linear substitutions. After 
providing some motivation, we start in Subsection 2.1 with an axiomatization of function 
algebras over a held K with a contravariant action of the matrix monoid (A4n(JL), •). These 
structures, subsequently named K-hierarehies (Dehnition 2.1), serve as the basic building 
blocks for the partial integral operator rings—the latter act on iP-hierarchies, and they 
take their coefficients from among them. Every iP-hierarchy induces a lattice of subalge¬ 
bras, which provide an algebraic formulation for functions depending on certain prescribed 
variables. In Subsection 2.2 we add a suitable collection of Rota-Baxter operators to a 
iP-hierarchy to arrive at the pivotal notion of a Rota-Baxter hierarchy (Dehnition 2.3). 
Its core ingredient is an algebraic formulation of those cases of the substitution rule of 
integration that are needed for obtaining normal forms in the operator ring. The basic 
properties of Rota-Baxter hierarchies are then developed (Lemmas 2.7, 2.8, 2.9). The 
central example of a Rota-Baxter hierarchy is the collection of all smooth multivariate 
functions (Example 2.5). 

The introduction of coefficients is studied in Section 3. Ignoring hrst the Rota-Baxter 
structure, we start in Subsection 3.1 by describing an algebraic structure that captures 
“separated” multivariate functions under linear substitutions. While the coproduct of 
the bialgebra provides the basic linear substitution x ^ x + y, proper linear combina¬ 
tions like X ^ X + \y can only be achieved via an additional “scaling” structure, thus 
leading to the notion of scaled bialgebra (Dehnition 3.1). Intuitively, the elements in 
such a scaled bialgebra are the univariate functions from which one builds all available 
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coefficient functions via linear substitutions. This is made precise by the induced hierar¬ 
chy (Dehnition 3.2), which is indeed a ii'-hierarchy (Proposition 3.3), and of course the 
smooth algebra—along with certain subalgebras—provides the prototypical model for this 
situation (Example 3.4). Next we combine scaled bialgebras with Rota-Baxter operators 
(Subsection 3.2). This requires a certain compatibility relation, which is in fact a special 
case of the substitution rule (Dehnition 3.5); the induced hierarchy of such a Rota-Baxter 
bialgebra is then a Rota-Baxter hierarchy (Theorem 3.7), and this holds in particular for 
the smooth algebra (Example 3.8). They form the appropriate coefficient algebras for the 
operator ring (Dehnition 3.9). 

Section 4 is devoted to the construction of the operator ring, focusing hrst on the iden- 
tihcation and verihcation of fundamental operator relations (Subsection 4.1). From the 
perspective of rewriting, the ehect of these relations can be described intuitively as follows: 
In the hrst step (Proposition 4.2), one normalizes one-dimensional integrators —meaning 
integral operators with subsequent linear substitution and coefficient—by ensuring that 
their associated linear substitution is “minimal” (for integration along the Xj-axis: identity 
matrix except beneath the i-th diagonal element). Then the normalized one-dimensional 
integrators are put in ascending order (Proposition 4.4), and repetitions are eliminated 
(Proposition 4.7). The monomials in the resulting normal form will therefore consist of a 
strictly ascending list of one-dimensional integrators, as will be shown later (Theorem 4.11). 
In the hnal Subsection 4.2 we dehne the desired ring of partial integral operators and linear 
substitutions (Dehnition 4.8). By its very dehntion, it is clear that it has a natural action 
on the given Rota-Baxter hierarchy (Proposition 4.9). The corresponding rewrite system 
is shown to be Noetherian (Theorem 4.10), meaning the computation of normal forms 
will always terminate. We conclude by conjecturing that the rewrite system is moreover 
confluent (Conjecture 4.12), meaning the normal forms are in fact unique. 

In the Conclusion we provide some pointers to possible future work. On the one hand, 
we propose some thoughts on the next obvious steps towards building up an algebraic and 
algorithmic theory of LPDE boundary problems. On the more speculative side, we discuss 
also some ideas about potential generalizations and extensions of our theory. 

1.3. Notation. The set N of natural numbers is assumed to include 0. If X is any set 
such that 0 G X, we write = {x & X \ x 0} for its nonzero elements. 

The vector space of m x n matrices over a held K is denoted by K'f, where m = 1 
or n = 1 is omitted. The unit vectors of iC" or Kn are written as e*; it will be clear 
from the context what n is, and whether rows or columns are intended. We denote the 
corresponding matrix rings by M.n{K) = Xj). The n x n identity matrix is denoted 
by In- We write Cij G Xin{K) for fhe matrix units, with the only nonzero entry 1 in 
row i and column j. Similarly, dj ^ = In + — l)^^ denotes the scaling matrix, which 

is equal to except that its i-th diagonal entry is A. The vertical composition of two 
matrices M E K^, N G is denoted by M © X G 

Given a K-algebra X, we write AlgHom^(X) for the monoid of JP-algebra endomor- 
phisms on X. The opposite of a monoid or algebra X is denoted by X*. If S' is a 
semigroup, we use the notation K[S] to refer to the semigroup algebra over K. In this 
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paper, all (co)algebras are assumed to be (co)commutative, so in particular bialgebras 
and Hopf algebras are both commutative and cocommutative. In contract, the word ring 
shall designate unital rings that are not necessarily commutative. We apply the usual 
notation [S,T] = ST — TS for the commutator of ring elements S and T. 

In this paper, all Rota-Baxter operators and derivations are of weight zero. Since we 
think of them as inspired from analysis, we write them as J and d. These symbols (along 
with their embellished variants) will always be used in operator notatation. For example, we 
write {dxf){^^g) rather than dx{f)^^{g)- Prodncts—usnally denoted by juxtaposition— 
have precedence over operators, so is to be parsed as and as ruTs)- 

Note that derivations like dx are indexed below to indicate their origin from in contrast, 
integrators like are indexed above as a reminder of 

In the sequel, we will use the standard variables Xi,X 2 ,X 3 ,... for dehning functions of 
arbitrarily (but finitely) many variables. So a definition like f{x 2 ) '■= X 2 should be carefully 
distinguished from f{xi) := Xi since the first denotes the a: 2 -projection {xi,X 2 , ■ ■ ■) ^ X 2 
bnt the second the Xi-projection (xi, X 2 ,...) i-a xi. Sometimes it will be convenient to use 
the abbreviations x = xi,y = X 2 ,z = X 3 . Since we do not use x as a shorthand for the 
seqnence (xi,X 2 ,X 3 ,...), this will create no confusion. In the scope of this paper we will 
only deal with linear snbstitntions, hence we will usnally drop the qualihcation “linear”. 

2. Rota-Baxter Hierarchies 

We start by bnilding np the basic operational domains—the algebras which the prospec¬ 
tive ring of partial integro-differential operators is to operate on. Such domains will be 
called Rota-Baxter hierarchies since they encode a notion of mnltivariable integration in 
conjnnction with an action of the corresponding substitntion monoid. 

2.1. Background to the Concept of Rota-Baxter Hierarchy. As a motivation to the 
general definition, let us first look at the classical setting of multivariate smooth functions. 
In the latter case, we are thinking of fnnctions /: K,"' —)■ R, for any arity n > 0, with the 
nullary ones denoting constants. To simplify the book-keeping we pass to the direct limit 

C°°(R°°) ;= IJ C'°°(R^) 

n>0 

of functions depending on arbitrarily (but finitely) many real variables. 

Thus we think of iF := C°°(R°°) as an algebra with the ascending filtration C Jm+i 
given by = C'°°(R"), so we can set up the action of arbitrary-sized real matrices. An 
ascending sequence of algebras {J^n) will be called an ascending algebra. 

For the rest of this section, let iF be a field of characteristic zero so that Q C K. We 
write M.{K) for the monoid of all cxo x cx) matrices M: N x N —)■ JF that can be written 
in the form M = I + M, where I is the 00 x 00 identity matrix and where M is any row 
and column hnite matrix (meaning all rows except finitely many and all columns except 
hnitely many are zero). We call Ad(iF) the (linear) substitution monoid over K. Note that 
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it has the natural ascending hltration 

M{K)=\J {l + Mn{K)), 

n>l 


where we use the embedding AiniK) Ain+iiK) that sends M to i). We will also 
identify hnite matrices M G Ain{K) with their embedding E M(K). In particular, 

we regard scalars X E K as Aen + 622 H— ■ G M(K) rather than XI. In the sequel, we will 
also need the descending chain 

M(K) = Mo(K) D Mi{K) D M 2 {K) D ■ • • , where 


Mn{K) 


0 

M 


M eM{K)^ C K[M] 


is the subring of matrices acting trivially on K"-. A rectangular matrix M E Kl with r < s 
is identihed with the corresponding square matrix M E obtained from M by adding the 
unit vectors e^+i,..., as additional rows. On the other hand, given M E Kl with r > s, 
we identify M with M E Kl by adding zero columns. In particular, any row v G K^ may 
be viewed as a matrix n© 62 © • • • ©e„ G it'”, and adjoining a column vector w E K'^ to the 
identity matrix yields the square matrix matrix {In,v) © e„+i G M.n+i{,K). As mentioned 
above, all square matrices are further embedded into K[M.] via Ain{K) Ain+i{K). 

In the classical setting A' = ]R, the canonical action of monoid Ai(lR) on the ascending 
R-algebra 0°°(]R°°) is dehned as follows. For given M E Aifc(R) and / G C°°(JR!‘), we 
set n = max(/c, 1) so that we may take M E Ad„(]R) and / G C'°°(]R”) via the corresponding 
embeddings. Then we think of M as effecting the change of variables 


Xi = miiXi H-h minXn, 


Xfi 7 TIyi\X\ © * * * © 

and we dehne the action AT„(R) x (^““(R”) -E- 0°°(R”) by {M,f) eA / o M. Via the 
embeddings, this yields the desired action W1(R) x C'°°(R°°) —)■ C°°(R°°). 

Let us now turn to the general case of a contravariant monoid action K[Ai] y. IF ^ IF, 
meaning a monoid homomorphism K[Ai]* -E AlgHom^(J^), where K[M]* denotes the 
opposite monoid of K[Ai'\. We require the following natural compatibility condition. For 
given M E M(K), we write M* for the induced mapping IF IF. Moreover, we write 
for the n-th cut-off substitution, whose hrst n rows coincide with those of M while the 
subsequent ones are e„+i,e„+ 2 ,.... In C'°°(R°°) this means substitutes only in the 
hrst n variables while leaving the remaining ones invariant. We call the action straight 
if M*f = for all M E M.{K) and / G IFn- The canonical action on C'°°(R°°) is of 

course straight while for example shifting the hltration to JV := C'°°(R”’''^) leads to an 
action that is not straight. 

Another crucial property of the classical example C'°°(R°°) is that evaluation of a func¬ 
tion / G C°°(R”) at Xn = ^ G R leaves a function in C°°(R”“^). In the general case we 
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shall require this only for evalation at .^ = 0, which can be described as the action E*, 
where En := / —G is the n-th evaluation matrix (at zero). Adding this require¬ 

ment to straightness, we arrive at the following axiomatization of multivariate functions. 

Definition 2.1. An ascending A'-algebra (Jm) is called a K-hierarchy if it has a straight 
contravariant monoid action of Ai{K) such that M*{En) C En for all M G Ain{K) 
and E*{En) C En-i- We write E for the direct limit of {En)- By abuse of language, we 
refer also to as a hierarchy. 

In detail, a contravariant action means each M*: is a homomorphism of K- 

algebras that restricts to a homomorphism M* : En whenever M G M.n{K). More¬ 

over, we have the usual action axioms I* = Ij- and {MM)* = M*M*. Note also that we 
have assumed Eo = K, which implies that the action on Eo is trivial (since it hxes 1 G K). 

In a A'-hierarchy {En), we can dehne the following dependency lattice for expressing 
that some functions depend only on certain variables Xa^, ■ ■ ■ ,Xak- For convenience, let us 
identify strictly increasing sequences ai < • • • < with hnite subsets {ai,..., C N; 

we will use the notation (ai,..., ak) for either of them. Furthermore, we shall iden¬ 
tify permutations tt E Sn with the permutation matrices (7 r(ei),..., 7r(eji)) G M.n{K). 
Let Soo := UnGN group of all permutations with finite support (those fixing all 

but finitely many elements of N). Then we have an embedding Soo ^ J^{K). However, 
note that the action on column vectors is accordingly 



where tt is the inverse of tt. We introduce the iF-subalgebras 

En.:={fEE\ Mf G Ek} 

where vr := tTq, : N —)■ N is the unique permutation with hnite support sending j to aj 
for j = 1,..., /c while restricting to an increasing map N \ {1,..., /c} —)■ N \ a. (In fact, the 
action of tt outside of {1,..., fc} is immaterial because of the straightness assumption.) By 
convention we set E% = E^. One checks immediately, using the straightness of the action, 
that a I—)■ Ed is a monotonic map (in the sense that it preserves inclusions). Hence we may 
extend it to arbitrary a C N by dehning 


OO 

‘Fq, = 

n=\ 

This yields a complete bounded lattice {Ea) with Ea^Ey = Ea<oy and EoEEy = Eany, with 
bottom element Eq = K and top element Am = E. Moreover, the lattice is complemented 
with E'd = FisjXa- Intuitively, Ea captures those functions that depend at most on the 
variables specihed in a and E^ those that do not depend on these variables. 
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The usual substitution notation /(MnXi + • • • + MinXn, • • •, M^iXi + • • • + MnnXn) can 
be viewed as a convenient shorthand for M*{f), for a given substitution matrix 


M 



and f E J^n- While we will not employ this notation in the present paper (a more suitable 
notation is introduced in Section 4.1), it is certainly useful in a computational context. 

The next step is to add Rota-Baxter operators and to describe their interaction with 
substitutions. This will lead to an algebraization of the well-known substitution rule for 
integrals. For basic dehnitions and properties of Rota-Baxter algebras we refer to [17, 18]. 

Let R be a ring containing Q, and let (R, P) be a Rota-Baxter algebra over R. Then 
we call {P, P) ordinary if P is injective and Im(P) -(- R = R as R-modules. This is an 
algebraic way of describing P as an integral operator on “univariate functions”. In fact, 
we get an ordinary integro-differential algebra {P,d,P), where d: R —?■ R is the unique 
derivation that sends R(/) -I- c to /, for arbitrary / G R and c E R. Hence 1 — R o d is 

an algebra homomorphism R —)■ R, which we call the associated evaluation of (R, d, R); it 
is the projector corresponding to the direct sum Im(R) -(- R = R. 

Having an ordinary Rota-Baxter algebra has the added beneht of having the polynomial 
ring at our disposal. This holds for all ordinary integro-differential algebras [ 8 , Prop. 3], 
but we give an independent proof here that does not make use of the derivation. 

Lemma 2 . 2 . Let (R, R) he an ordinary Rota-Baxter algebra over R. Then x i-A R(l) 
defines an embedding (R[a:], ^ (R,-P) of Rota-Baxter algebras. 

Proof. Since (R[x], J*) is the initial object in the category of Rota-Baxter R-algebras [17, 
Cor. 3.2.4], there is a unique Rota-Baxter morphism t: R[a;] R, which clearly satis- 
hes i{x) = i'(Jgl) = R(l)- It remains to check that i is injective, so we show that i{p) = 0 
implies p = 0 for all polynomials p E R[x]. We use induction on the degree of p. The 
induction base p G R is trivial since by dehnition l acts as the identity on R. Now assume 
the claim for all polynomials of degree less than k > 0, and take p = po + Pix pkX^ 

with l{p) = 0. Using the property R(l)* = i\ R*(l) from [17, Thm. 3.3.1], we obtain 

k k /k-i \ 

-Po = ^ + l)!pm^*(l) e MR), 

i=l i=l \ i=Q ) 

and Im(R) -(- R = R implies that po as well as the expression on the left-hand side above 
vanish. Since P is injective, this implies 

( A:—1 \ k—1 k—\ 

+ l)pi+iR J = + l)pi+iR(l)* = ^{i + ly.pi+iPfil) = 0, 

1=0 } i=0 i=0 

and the induction hypothesis yields (i-l-l) Pi = 0 for alH = 0 ,..., k—1, and hence p = 0 . □ 
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For an ascending algebra (Jm), it is natural to require an infinite collection of Rota-Baxter 
operators that we shall write as Since we think of IFi as univariate functions, we shall 
require that is an ordinary Rota-Baxter operator over K = J^q. Analogous assumptions 
are imposed for so that is an ordinary Rota-Baxter algebra over IFn-i- For the Rota- 
Baxter operators we shall now postulate the substitution rule for integration, which 
we shall only need for certain particular substitutions. 

We introduce the following two special matrices. For i < nwe define the general transvec- 
tion (= horizontal shear) in the Xj-direction as 


/I 


\ 


( 1 ) 




Vl 


1 


Vi-i 1 Vi+I 
1 


Vn 


^ row i, 


\ 


1 / 


for a vector conveniently written as n = (ui,..., Vi^i, Uj+i,..., VnY € Similarly, we 

define the eliminant (= subdiagonal vertical shear) 


( 2 ) 


/I 


Li{w) = 


1 

Wi+i 1 


\ 


^ row i + 1, 


V Wn 1 / 


for a vector written as w = (tCi+i,... ,tCn)^ G K^~\ When using the abbreviated vari¬ 
ables (x, y) = (xi, X 2 ), we shall also write T^^v) and Ty{v) for Ti(n) and T 2 (n), respectively, 
and similar abbreviations will be in force for the eliminants. We note also the composition 
rule Li{w) Li{w) = Li{w-\-w) so that Lf^{v) = Li{—v) for eliminants in the same direction. 
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However, for distinct directions we have the rules 


/I 


\ 


1 


1 

Wi+i 1 


<(— row i + 1 


(3) Li{w)Lj{u) 


1 

Uj _|_ 1 1 


row j + 1 


V 


(4) Lj{u) Li{w) = Li{w') Lj{u) 


with i < i and w' := Lj_i(u)w G K"‘~\ In particular, we note that the Li{w) and Lj{u) 
do not commute. 

We can now introduce the key structure that we shall use as an algebraic model of 
multivariate functions under integration and substitution. Roughly speaking, we shall need 
the subsitution rule for integration just for the following three matrix types: scalings, 
transvections, and eliminants. 

Definition 2.3. A Rota-Baxter hierarchy consists of a A'-hierarchy and 

commuting Rota-Baxter operators (J^") that satisfy the following axioms: 

(a) We have ^ and for n < m. 

(b) Every J'") is an ordinary Rota-Baxter algebra over iFn-i with evaluation Ef. 

(c) We have r* for the transposition r = (i j). 

(d) For A G and v = (0,n') G iF © we require the substitution rules 



(5) 

( 6 ) 
(7) 


We call (5), (6), (7) the diagonal, horizontal and vertical substitution rule, respectively. 

Note that (5) describes the effect of scaling, (6) is essentially a one-dimensional substi¬ 
tution rule while (7) is the signihcant part of the two-dimensional substitution rule, which 
will be used for collapsing multiple integrals along the same direction. It turns out that 
all other instances of the substitution rule needed here can be inferred from the three 
instances (5)-(7). In particular, we note immediately that they imply the slightly more 
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general cases of the three substitution rules 


( 8 ) 

( 9 ) 

( 10 ) 


r^i 


rufe)’ = (i-u)u(oT 

S^‘Li{e, + vrr = L-_;,{vr[L,(e,r,p]r*‘r+i(v'y 


for any i,j > 0, as one may readily check using conjugation by transpositions. In the 
vertical rule (10) we require v = (0,n') G 0 and j < i. 


2.2. Examples and Properties of Rota-Baxter Hierarchies. The vertical substitu¬ 
tion rule (7) can be formulated in the following equivalent way, which is more symmetric 
and perhaps more natural (but less economical for our present purposes). This can be use¬ 
ful for proving that something is a Rota-Baxter hierarchy (as we will do in Theorem 3.7). 


Lemma 2.4. Let A C be an arbitrary finite index set with minimal element A G A and 
complement A' = A \ {A}. Then 

(11) jn.( 5: e.) ’r = ioli (E n (E ^<-a) ■ 

ieA ieA' ieA' 

is equivalent to (7), assuming all other axioms of a Rota-Baxter hierarchy (JTi, r )nGlN • 


Proof. It is easy to see that (11) is necessary. Indeed, let n be the largest index of the 
set A. We may assume that A = 1 since otherwise we can set up a permutation vr G S'„ such 
that 7r(A) = 1; then conjugation of (11) by tt* will ensure this. Now setting v' = Yhi^A' 
in (7) immediately yields formula (11). Note that the appearance of the ej_A = ej_i in the 
latter formula corresponds to the transition from v G 77”“^ to v' G 77”“^ in the former. 

For proving that (11) is sufficient, let v = (0, n') G 77 © 77"'“^ be given, and let A' C 
denote the set of all those i with v[ 7 ^ 0. Furthermore, set A = 1 and A = {1} U A'. We 
use the scaling matrix S = S 5 ■ ■ ■ Sn & j^nxn 


^*+2 = 


di+2,vr if Vi 7^ 0, 

In otherwise. 


for 1 < i < n — 2. Then we have Lx{ei v) = S Lx{J2ieA ^ ^ hence 

+«,)• = (©)•... (SU)-LhE' 


iGA 


Q* Q* 

*^3 * ‘ 


Since the S*, (S'- ^)* G A42{K)* commute with , substitution into (11) yields 


rn(£!.+«)• J'=(s-‘)* Lp (E e.-i) ■ [LMT, n s’lJ e 


iGA' 


, V-i—i 


iGA' 


s* 


{S-'r Lp ( E e.-i) A- [LMT, f] f (S-fU ( E ei-i) 's’, 


ieA' 


ieA' 
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where in the last step = 1 was inserted after and S* was pushed left because 

the S* also commute with and L^^ei)*. Then we obtain (7) since clearly S = 
L~^[v') as well as SLy{'^) S~^ = Ly{y'). □ 

Let us now look at the most important example of a Rota-Baxter hierarchy—the algebra 
of smooth functions in several variables. This hierarchy contains also several important 
subhierarchies, in particular the analytic functions. 

Example 2.5. In the classical setting we have iFn = C'°°(]R”) with Rota-Baxter operators 

rxn 

Jo 

which clearly satisfy the mapping properties required in Item (a) above. Moreover, it is 
clear that every is injective, and we have Fn = Fn-i + Im(J'^") since 

f{xi, ...,Xn)= fixi . . . , Xn-l, 0) . . . , Xn- 1 , 0 df, 

and since 0 = c(xi,..., Xn-i) + /q^" /(xi,..., Xn-i,0 df implies / = 0 upon differentiating 
with respect to Xn- Thus every {Fn, is an ordinary Rota-Baxter algebra over Fn-i with 
evaluation Ef:Xn^ 0, as required by Item (b). The transposition property of Item (c) is 
clear, while (8) follows by the substitution .^ = A.^ in the integral X 2 , ■ ■ ■) df. 

Now let us consider (6). Writing T = F^ei), we have 


fXl 


f{f + X2,X2,X3,...)df = 


rxi+x2 




fX2 


(*X2 


X2, X3,...)df={l- E*) T* X2, X3,...) 


j''^T*f{xi,X2,X3,. . = 

pXl+X2 

= / f{^F2F3,---] 

Jo Jo 

where the second equality employs the change of variables = ■C + 2 ^ 2 - 

Finally, let us verify (7). Using the abbreviation Zk...i = Zk, ■ ■ ■, zi {k < 1) for extracting 
and manipulating partial vectors with the obvious meaning, we have 

L 2 {v’yp^ L^ei + v)*j''^f{xi, ...,Xn,...) 


= L2{vy 


rxi 


/ /(C, X2 + T], X3,„n + V3,„nri, . . .) df df] 

) Jo 

fXl l'Xl+X2 


= L2wr 

Jo 

nxi l•x\+X2 


fif, X3,„n + V3,„n{r] - X 2 ), ■ ■ ■ ) dt] df 


C+X2 


10 


>0 


^+X2 

xi rxi+X2 
<0 


fif, V, X3...n + V3,„nV, ■ ■ ■ ) df] df 

m r^+x2 


... dr] df 


... dr] df 


where in the second equality the integration sweeps are swapped and the substitution t] = 
X 2 + r] is applied. It is easy to see that the hrst summand is ^ 2 ,( 61 ) ^rriyi v'Yf and 
the second is Lx{ei)*Ly{vy f. 
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This concludes the proof that C'°°(]R°°) is a Rota-Baxter hierarchy over ]R. A similar 
proof will also work for the analytic subhierarchy C'‘^(1R°°). This may be extended to 
complex variables as follows. Let J be the path integral from 0 to x G C. If / G C‘^(C°°) 
is a multivariate holomorphic function, we use the conjugates = (1 (1 i)* for 

creating the hierarchy of Rota-Baxter operators. One sees immediately that is 

a Rota-Baxter hierarchy over C. (There is also an intermediate case where one considers 
only comp lex-valued functions on R", but allowing complex substitutions. For example, 
e“ = cos(x) -|- i sin(x) may be interpreted in that way.) 

Remark 2.6. The classical example provides a convenient geometrical interpretation of 
the three substitution rules (5), ( 6 ) and (7). The diagonal rule (5) describes the natural 
contravariant behavior when stretching or shrinking a coordinate axis. 

The horizontal rule ( 6 ) says that integrating over the segment [0,x] x {y} after hor¬ 
izontally shearing it to [y, x + y] x {y} can be achieved by integrating over the whole 
sweep [0,x + y] x {y} minus the surplus [0,i/] x {y}. Its effect is that all axis-parallel line 
integrals may be started from the axis, generalizing the usual rule . .dx = jo---dx- 

... dx on the x-axis to all parallel segments. 

The vertical rule (7) for n = 0 may be seen to decompose an integral over the triangle 
( 0 , y), ( 0 , x + y),{x,x + y) as an integral over the rectangle ( 0 , 0 ), (x, 0 ), (x, x + y),{0,x + y) 
minus an integral over the complement, namely the trapezoid ( 0 , 0 ), (x, 0 ), {x,x + y), ( 0 , y). 
Parametrizing the line segment from (0,i/) to (x, x + y) as rj = s{^), the former integral 
is JoJ^• . dr] d^, the latter is Jq ... dr] d^, so the effect is again to offset all integrals 
from the axes. The case n 7 ^ 0 is similar except that the triangle is now tilted against the 
xi/-plane, but the same decomposition as before appears in the orthogonal projection. 

Finally, we should point out that the special case A = —1 of (5) means that (—1)* and 
anti-commute. Thinking of Dehnition 2.3 as an axiomatization of the substitution rule of 
integration for the case of linear substitutions A G ATn(R), the proper generalization to 
general “spaces” must be the signed integral over oriented manifolds rather than the so- 
called area formula of measure theory [14, §3.3] [42, Thm. 2.6]. The latter would introduce 
absolute values whose algebraic treatment would be considerably more awkward. 

Clearly, the notion of Rota-Baxter hierarchy gives rise to a category, with the morphisms 
respecting the linear action as well as the Rota-Baxter operators. In connection with this, 
several other notions can be built up in a natural fashion but in the present context we 
shall not need them. Let us only mention that a Rota-Baxter subhierarchy means that 
corresponding hltered components are Rota-Baxter subalgebras in the usual sense and 
that the linear action of iF[A4] coincides on them. 

We proceed now by stating a few elementary conseguences of the axioms. Though most 
of these are fairly obvious in the classical setting, we have to prove them here on the basis 
of our axioms for making sure that these include all the essential requirements for the 
theory to be developed. 

Lemma 2.7. Every Rota-Baxter hierarchy (JTi, satisfies the following properties. 
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(a) For any a = (ui,..., there is an embedding 

Xaj Xaj ■= J “^1, 

and we have 7r*p(a:ai,..., Xq^,) = ..., x-^i^ak)) all permutations tt o/(ai,..., a^). 

For 71 E Sn and i < n we have 7i* r = r '^71* with j := 7i{i). In particular, 
all IF{i) are conjugates of and hence ordinary Rota-Baxter operators. 

(c) We have f for all c E and f E T. In particular, = cx„. 

(d) The embedding ia of Item (a) is a homomorphism of Rota-Baxter algebras in the sense 

that La o o La for j = 1,..., k. 

(e) If M E K[M] vanishes in the i-th column, then M*{IF) C JF. 

(f) We have E* = 0 for all i > 0. 


Proof, (a) It suffices to consider a = (1,... ,n) since restricting to a subset of a induces 
another embedding, and permutations may be restricted accordingly (fixing the elements 
outside the subset). 

By dehnition, := La is a homomorphism of J’T-algebras. We show that it is injective 
by induction on n, with the base case n = 1 covered by Lemma 2.2. Since Ea = En is 
an ordinary Rota-Baxter algebra over En-i, applying Lemma 2.2 again yields an embed¬ 
ding l: (J-n-i[^n], ^ defined by i-A Xn. By the induction hypothesis 

we also have an embedding Ln-i'. K[Xi,... ,X„_i] ^ En-i defined by W i-A Xj (i < n). 
By the universal property of polynomials [5, Lem. 2.15], we obtain a unique JL-linear 
map Ln'. K[Xi,..., Xn-i\[Xrf\ -E En-, wlfich acts on coefficients via and sends X^ 
to Xn. By [5, Lem. 2.16] and the injectivity of tn-i, the map Ln is injective iff 


k 

(12) '^pii^xi,... ,Xn-i) Xn = 0 implies po,...,pk = 0. 

i=0 


For showing (12) note that the left-hand side is the image of ^n-i{Pi) under the 
embedding l, so that po, ... ,pk = 0 follows from the injectivity of l and in-i- 

Now let TT G Fn be an arbitrary permutation. Then Item (b) of this Lemma (whose proof 
below is independent) yields 


7l*p{xi,. .. ,Xn) =p(7r*J'''l,.. . ,7r*J"'"l) =p(J"'"(^^l,..., =p{ 


X, 


( 1 ), 


, X. 


r(n)) 


since clearly vr*! = 1. 

(b) The general conjugate relation follows from Item (c) Definition 2.3 since every tt E Sn 
is a product of transpositions. The last claim follows since = T* with T = (li) 
is clearly injective and for given / G E^i) we have t* f = fo-\- fi for some fo E K 
and /i G E(i), which implies / = /o + for fi = x*fi E E(i), and the decomposition is 
unique since the one for ^”( 1 ) = Ei is. 

(c) For a sufficiently large k > n we have f E Ek and c G E[a) with a := {1, ... ,k) \ (n). 
So if TT is any permutation taking to i, we have 7t*c G Ek-i. We choose a permutation 
with 7i{n) = k so that vr*/ G Ek. Then Item (b) of Definition 2.3 yields J'^'‘(7r*c • tt*/) = 
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(7r*c) If TT is the inverse of vr, left multipication by n* together with Item (b) of 

this Lemma gives the required identity. 

(d) Again it suffices to consider a = (1,..., n), so we show that o = r o for 
hxed n. In the case i < n we use hnite induction on i. The induction base i = 1 is then 
covered by Lemma 2.2 and Item (c) above, so assume i > 1. Then we have 


ir,. O 


rXi / k 
j=0 

rXi 


rXi 


,.o •'» 


n—ly 


j=0 

f* 'T' ■ k 


Pj{Xi,...,Xn-l)] =^xi / Pj{Xu...,Xn-l) 


/ Xi ^ rXi 

'^Pj{Xi,...,Xn-l)xi^= Oi„ 

j=0 


j=0 
k 


j=0 


where we have used the induction hypothesis in the third and Item (c) of this Lemma 
in the fourth equality. It remains to prove in o I„ = f” on = 

K[Xi,... ^Xn_]\[Xj^. To this end, recall that the embedding 1 ^ 
from the above proof of Item (a) is a Rota-Baxter homomorphism by Lemma 2.2. Hence 
we obtain 


pX„ / k 


Ln O 


j=0 


j=0 


niXi, ...,Xn-i)xi)=i„(J2 Pi(Xu ..., w_i) x’PKj + 1 ) 

k 

^Pj{xu...,Xn-l)Xi 


rXn / k 


^Pj(xi,...,a:„_i)xf,+V(i + 1) = k 

j=0 

k 


j=0 


O L 


^Pj{Xi,...,Xn-l)X' 
j=0 


3 1 - 

n 


OLr 


j=0 


where the Rota-Baxter property of l has been employed in the fourth equality. 

(e) We may assume M G A4n(X) and f ^ Xn for a sufficiently large n >i. Setting r = 
(i n) G Sn, we see that Mr has the last column zero so that Mr = {Mr) E^. Hence 
we have also T*M*f = [Mr)* f G by the evaluation property of hierarchies in 

Dehnition 2.1. Using the cyclic permutation 


7 = (zi—— It—)- z)g Rn-i, 

we see that ry is a permutation that sends fc to cifc for a = (1,..., i — 1, i -1- 1,..., n), so by 
the dehnition of the dependency lattice we may infer M*f G Xq, ^ from {T'y)*M*f = 
j*T*M*f G Xn-i- But the latter follows immediately from T*M*f G Xn-i since 7 G Sn-i- 

(f) Take f e Xn- li n < i then f E Xi and hence E*f = 0 follows from Item (b) of 
Dehnition 2.3. Otherwise let r be the transposition (n i) so that = 0 by the 
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same token. Composing this on the left by t* gives f = 0 by Item (b) of this Lemma 
since Ei = tE^t. □ 


The horizontal substitution rule (9) can be generalized to the following result about 
arbitrary transvections. 

Lemma 2.8. Let (Jti, be a Rota-Baxter hierarchy over a field K. If T = Tfiv) 

is any transvection along the Xi-axis with {vi,... ,Vi-i,Vi+i,... ,Vn) € Kn-i, then we 
have p T* = (1- E*) T*p. 


Proof. The general result follows from the case i = 1 by conjugation with (1 i) G Sn and 
Item (c) of Dehnition 2.3. Hence we may assume 7 = 1. If denotes the transvection (1) 
corresponding to v = {v 2 , , Vn) G K^-i, we have the obvious relation Ty = • • • Tv„en 

for the composition, so it suffices to consider := T„ for v = acj (j > 1). We may 
furthermore assume that a 7 ^ 0 since otherwise Ta = In, and the result follows by Item (f) 
in Lemma 2.7. If S' = dj^a denotes the corresponding scaling in Xj-direction with inverse S = 
dj,i/a, we have = STiS and hence 


p Tf = pS*TfS* = S*pTfS* = S*{1- EDppS* 

= (1 - u) t: r 


s’{i- 


where the second and fourth equality follows from Item (a) of Dehnition 2.3, the third 
from (9), and the last from EiS = SEi. □ 


The vertical substitution rule can also be generalized in a way that will become important 
in the next secion—allowing a coefficient function within the integral. For handling this 
kind of situation we use the technique of slack variables. Before going through the formal 
proof, it may be helpful to recall this technique from Analysis. In fact, the verihcation of 
the vertical substitution rule in Example 2.5 can be generalized as follows: 


L 2 {v')* p g{xi) Lpi + vypf{xi, ...,Xn) 

rxi pr) 

= L 2 {v')* g{r]) f{i,X2+V,X3...n + V3,„nV)dfdT] 

Jo Jo 

/•Xi pxi+X2 

= L 2 {vy / g{fi - X 2 ) f{^,f],X3,„n + V3,„n{fj - X 2 ))df]df 

Jo j (,+X2 

pxi pXl+X2 

= / g{xi,r],X3,... ,Xn,X2) f{^,r],X3,„n + V3,„nr])dr]df 

Jo J ^+X2 


|•X1+X2 rx\ I‘i+X2 

I ... dr] df — / ... dr] df 

0 Jo Jo 

The auxiliary function g G En+i is dehned as g{xi,X 2 , ■ ■ ■ ,Xn,Xn+i) ■= g{x 2 — Xn+i). 
Since the substitution induced by J := /„ © 62 acts as Xn+i •—t X 2 , it is now easy 
to check that the hrst summand is given by J*Lpi)*pg Ly[v')* f and the second 
by rpLpiYpgLyP) *f. Thus we recover the same rule as before, apart from the 
presence of the multipliers g,g and the back substitution effected by J*. The whole point 
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is that the slack substitution X 2 Xn+i allows us to temporarily “freeze” the variable X 2 
so that it is not affected by integration. However, there is a price to pay for this: The 
identity is now restricted to J^n rather than being valid over all . This is reflected by the 
change from /(xi,... ,x„,...) in Example 2.5 to /(xi,... ,Xn) in the above verihcation. 
Practically speaking, this is not an essential restriction: The slack variable Xn+i must be 
chosen large enough to prevent any conflict with the substitutions or integrals. Let us now 
prove the same result for general Rota-Baxter hierarchies. 


Lemma 2.9. Let J^"')neN be a Rota-Baxter hierarchy over a field K, and let g ^ 
be any coefficient function. Then on Tn we have the identity 

(13) f'sijlei + t.)T = ®ei+i)’[U(ei)'.r]r*‘9n+i(K')- 


for any 0 < j < i and v = (0, v') G iP* © iP"' * Here g := (cj+i — en+i)*g G Hn+i is 
the transform of the coefficient function. 


Proof. As before it will suffice to prove the case i = j = 1 since the general case can 
then be recovered via conjugation by transpositions. Let u G Hn be arbitrary but hxed. 
Since g G J^i, it is invariant under Lx(ei +x)* by the straightness of the action. Moroever, 
setting g := {—en+iYg G H(n+i) we have g = e„)*^, by straightness again since 

the hrst row of the matrix ( — e^+i © (0 /„)) Lx{—en) G is ei. Using this 

factorization we obtain 

^""g L^ei + = ^LYei + /""u 

= rK+i Lx{ei +v- CnYa Fm, 

where the hrst equality uses the multiplicativity of substitutions and the fact that G Hn 
is invariant under U*.,.! and Lx{—enY by straightness once again; the second equality 
follows from [Lx(ei + v),En+i] = 0 and (3), where L^^ei + n) G is embedded 

into via the hltration of Ai{K). Since g G H(n+i), we may move it in and 

out of the inner by Item (c) of Lemma 2.7, while En+i G A4n(E) and Item (a) of 
Dehnition 2.3 allows us to extract Ef_^^ from the outer f^. Then we apply (7) to obtain 

f'g L^(ei + = E*_f_^ + n - en)*J''gu 

= Ef^^ Lf^iv' - eY* [LxYYY r] r^yiv' - CnY 9U, 

and we observe that L~^{v' — eY) En+i = (/„ © 62 ) and Ly{v' — CnYd^ = g Ly{v'Yu, 

which may be verihed by a short calculation (again using straightness for the latter). □ 

The technique of slack variables gives considerable power to the notion of Rota-Baxter 
hierarchy. For example, choosing T = C'°°(]R"), the convolution of univariate func¬ 

tions -k: El X El ^ El can be dehned by 

/ A ^ := (Ji © eY*r{e^ - eyYf e*yg, 

which means {f x g){x) = j^f{x — y)g{y)dy, recovering the classical dehnition of the 
Duhamel convolution. Adding evaluations as in the ordinary case [35, §3], and limits as 
“evaluations at inhnity” [1, §4], one can dehne and study integral transforms (e.g. Fourier, 
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Laplace) in this algebraic framework. We will not pursne these topic in the present paper. 
For us, the main use of slack variables is to carry around coefficient functions within integral 
operators. This is a topic that we shall now investigate in some detail. 

3. Rota-Baxter Bialgebras 

Before we can build up the operator rings for Rota-Baxter hierarchies, we must hrst 
address the question of suitable coefficient domains. It turns out that such domains not 
only have an algebra structure (needed for composing multiplication operators) but also a 
coalgebra structure (needed for expressing basic linear substitutions) and a scaling action 
(in conjunction with the coalgebra this yields all linear substitutions). 

3.1. Conceptualizing Linear Substitutions via Scaled Bialgebras. Again it is help¬ 
ful to hrst look at the classical example C'°°(]R°°). Intuitively, we would like to simplify an 
integral operator like J^/(x, ?/), acting as u{x,y) ha f{x,y)u{x,y), by pulling out of the 
integral those parts of f{x,y) that depend only on x. For example, if f{x,y) = {x yY 
we would simplify 

(14) /^(x -F yf u{x, y) = u{x, y) + 2x ^y u{x, y) + JV u{x, y). 

However, this kind of simplihcation is not possible for integrands like f{x,y) = e^^ except 
if we are willing to use inhnite expansion like = ^^(l//c!) In the terminology 

of integral equations, f{x,y) = {x + yY is called a separated^ kernel, as opposed to the 
non-separated kernel f{x,y) = e^^. Since we would like to refrain from using inhnite sums 
(and hence topology), we will only allow separated kernels as coefficients in this paper. 
This will be made precise in Dehnition 3.9 below. 

The expansion step (14) can be understood as the substitution x x + y on g{x) = 
to yield f{x,y) = x"^ + 2xy -|- y"^. In other words, we have used the coproduct A: K[x\ —)■ 
K[x\ ® K[x\ = K[x,y] dehned by A(x) = x(8)l-|-l{8)xso that / = ^{g). Note that A 
is an algebra homomorphism; in fact, K[x\ has the structure of a bialgebra. For general 
properties of bialgebras we refer to [9, 12, 41] and to [17, §2]. As we will make precise later 
(Dehnition 3.5), the coproduct interacts nicely with the Rota-Baxter structure. 

Before we look at this interaction in more detail, it is apposite to focus hrst on the substi¬ 
tution structure of those bialgebras that provide “separated kernels” like the paradigmatic 
example K[x, y] above. It turns out that we can build a iF-hierarchy from a given bialge¬ 
bra "H from just one more ingredient, which we call scaling-, an action of the ground held K 
that will be extended to an action of the full matrix monoid K[M] on the tensor algebra 
over R. In the classical example, this is the action /(x) ha /(Ax) for a function / G C°°(]R) 
and a scalar A G IR^. 

We formulate the basic properties of such an action in terms of the convolution product 
that we denote here by 4-. As stated in the Introduction, all bialgebras are assumed to be 
commutative and cocommutative. Recall [17, Thm. 2.3.4] that for a iF-bialgebra "H with 
product V, unit 1, coproduct A and counit I one dehnes the associative and commutative 

^In Analysis, the term degenerated is most commonly used. From the viewpoint of Algebra, however, 
this term sounds too drastic so that we prefer the expression separated. 
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operation 4- on vector space endomorphism by setting f 4- g = V {f ^ g) A. If f and g are 
bialgebra endomorphisms, then both f 4- g and fog are. Writing BialgHom^('H) for the 
set of bialgebra ii'-endomorphisms, this yields two operations 

4-, o: BialgHom;^(?f) x BialgHom^('H) ^ BialgHom^(?f), 

which are clearly also ii'-linear. Moreover, one checks that {f 4- g) o h = [f oh) 4-[g oh) 
and h o (/ 4- g) = (ho/) 4 - (h o ^i). Hence (BialgHom^('H), 4 -, o) is a commutative 
unital semiring. The neutral element with respect to 4 - is the composite e ;= 1 o J of the 
unit 1: if' —?■ "H and the counit 1 ■. hi ^ K. Of course, id-^: if —)■ if is the neutral element 
with respect to o. 

Note that a (multiplicative) group action x "H —> "H is by dehnition a monoid 
homomorphism {K^,-,l) ^ (AlgHom^(if), o, id), which we shall write A i-A A*. Clearly, 
we may extend it to a monoid homomorphism K —)■ AlgHom^(if) by setting 0* := e. If 
the A* are bialgebra homomorphisms, we refer to x hi ^ BialgHom;^('H) as a group 
biaction. For a scaling, we want this homomorphism to respect the convolution product. 

Definition 3.1. Let ("H, V, 1, A, I) be a JF-bialgebra. Then a group biaction x "H —"H 
is called a scaling if the map (iC,+, •) (BialgHom^('H), 4-, o), given by A t-A A* is a 
semiring homomorphism. In this case, we call "H a scaled bialgebra. 

In the sequel, we shall suppress the notation o for the composition of endomorphisms. 
From the dehnition, we have the additive law (Ai + A 2 )* = A^ 4- A 2 = V(A* <8 A^jA. By 
induction, this generalizes immediately to 

(15) (Ai + ■ ■ ■ + A„)* = V" (A* (8) • • • ® a:) A" 

for all n > 0, where V" and A"' denote the obvious iterations (see below for the dehnition). 
Observe that the image K* C BialgHom^(?f) of the scaling homomorphism A i-A A* is 
automatically a field. Hence we have a held homomorphism K —)■ iC*, and the group 
action x is faithful. Furthermore, note that is a Hopf algebra with the 

antipode (—1 )*'. T-i T-L since this is clearly the convolution inverse of id^ = 1*. 

We write SfifH) = ®i>QHi for the tensor algebra with grades I-Lq := K and "Hj = "H®* for 
f > 0. The product on "H is denoted by V: "H (8) "H —)■ "H, its iterations by Hn -A H- 
Likewise, we write A": "H —)■ 'Hn for the iterated coproduct dehned by A^ = id^ and 
^n+i _ (gj Iterating coassociativity one obtains 

(16) (A™1(8)---(8)A”^")A^ = 

for all rui,... ,mk > 0 as in [28, Lem. 1.1.9(2)] with slightly diherent notation. Set¬ 
ting mi = 2, m 2 = ■ ■ ■ = mn = 1 yields the alternative recursion = (A (8) id”~^) A” 

for computing the iterated coproduct. The counit axioms imply 

(17) (e®(*-i) 0 id (8) E®(”-*)) A" = 0 id (8) 

where the right-hand side is the embedding "H —> Hn dehned by / i-A l®h-i) 0/0 
One checks also the iterated scaling commutations A*V" = V"'(A*(8)- • -OA*) for the product 
and A^-A* = (A* <8) • • • (8) A*)A"' for the coproduct. 
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The induced maps on the n-fold tensor products are then V"'®” = (V"')®"': T^®” -A l-iri 
and A"'®"' = (A")®"': T-in —t "H®”. Note that is a bialgebra with (iterated) product 
V"'®"'r and (iterated) coproduct rA”®", where r := r„: "H®” "H®" is the transposition 

(/ll (8) • • • <8 /in) (8) • • • <8) (/nl (2) • • • O fnn) ^ {fll ® ® fnl) ®---® (/in ®---® fnn)- 

If a = (oi,..., On) is a column vector in K'^ or a row vector in Kn we write a® = • ■®an 

for the corresponding scaling map on "Hn- Likewise, for a matrix A G M.n{K) we write 
A® = Af, ® ■ ■ ■ ® A®, for the scaling map on "H®”. If A denotes the transpose matrix, we 
have rA® = A®r or, in other words, r(Af, 0 • • • (8) A®,) = (Af^ 0 • • • 0 ^fn)'^- 

We can now dehne the matrix action in terms of the coproduct and the scaling action. 
Note the appearance of a single r in the expression for M* below: It allows us to interpret 
the left expression V”®”r as the product map on Tin or, via rM® = Mr, the right 
expression tA"®” on T-in as the corresponding coproduct map. 

Definition 3.2. Let ("H, V, 1, A, I) be a scaled bialgebra. Then the induced matrix action 
Mn{K) "Hn, Written (M,/) ^ M*f, is dehnded by M* = 

algebra io('H) together with the corresponding action K[AA] x S){l-L) —?■ is called the 

induced hierarchy for "H. 

If we dehne the action a*: H ^ Tin of a row a = (oi,..., a^) G by cl* = a® A”, the 
hnite matrix action A4n(K) x Hn —> 'Hn can be written as 

(18) M*{h 0 ... 0 /„) = (Mlfi) ■ ■ • (MUn), 

where juxtaposition on the right-hand side denotes the iterated product in The ex¬ 
tension of the hnite matrix action to K[M] x ij('H) —)■ io('H) is to be understood as in 
Section 2: Every matrix of K[AA] is of the block form (^?) with A G M.r{K) and / 
the oo X oo identity matrix, while every / G is a hnite sum of tensors of the 

form /i 0 • • • 0 /s 0 1 0 1 0 • • •, so it suffices to choose n as the maximum of those r 
and s and then use the action AT„(iL) x d-Ln -t- 'Hn- Let us now make sure that the 
induced hierarchy deserves its name. 

Proposition 3.3. Let (H, A, J) he a scaled bialgebra. Then the induced hierarchy S){H) is 
a K-hierarchy in the sense of Definition 2.1. 

Proof. The tensor algebra S)(H) is clearly an ascending algebra (Hn) with Hn = "H®” and 
direct limit Sj(H). Moreover, it is clear from the dehnition of the induced matrix action 
that M*(Hn) ® Hn for all M G Mn(K). 

Let us show that Ef(Hn) ® Hn-i- We write Lj := e 0 • • • 0id^ 0 • • • 0 e; Hn —t Hn with 
all entries E except for the j-th, which is id^. Likewise, we write rj := e®- ■ ■®e: Hn -A Hn 
for the tensor map having only E. Then E® = li® ■ ■ ■ ® Ln-i 0 p, and (18) yields 

0 • • • 0 /n-l 0 fn) = (mA"/i) • • • (i„_iA"/n-l)(r7 A"/„). 

Since e: H ^ K C H, we have rj AP fn G iL®” = K, so it suffices to show CjAP fj G Hn-i for 
all i < n. Writing lj : Hn-i -A Hn-i for the corresponding tensor maps with e everywhere 
except for the j-th entry id^, we have Lj = lj ® e and therefore 

LjA"' = (lj 0 e)(A”“^ 0 id^)A = (ijA"“^ 0 id^)(id-H 0 e)A = 0 1, 
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where the last identity uses the dehning property of the counit in "H. But this clearly 
implies that G 'Hn-i for all j < n, as was required. 

Next we prove that M.n{K) x Tin 'Hn is a contravariant monoid action. For any 
matrices M, M G Mn(-A) and any function / := /i ® • • • <8 /n G "Hn we must show 

(MM)*(/i ® ® /„) = «) • • • (8) /„). 

We start by computing M*(/i < 8 ) • • • < 8 ) fn)- Using sumless Sweedler notation, we have 

A”®"/ = (/l,(l) 0 • • • (8) /l,(n)) 0 • • • 0 (/n,(l) 0 • • • 0 fn,{n)) = <S> <S> /i,(jF 

* j 

where in the last expression (as in the rest of this proof) all indices range over { 1 ,..., n}. 
We have then M® A"'®"'/ = (^- (^j M*^ and further 

(19) V7 = (SinMl/,,|,|, 

j 1 

with denoting the iterated product of Ti. Using the fact that A” is a bialgebra morphism, 
this implies A"'®"'M*/ = Now observe that the M^,..., M®, are 

algebra morphisms so that 

A)r®A’'®”M7 = (g) n (® (A7i,A = n ® (® (® /..OjA ■ 

= n f0 0 {AUjMjk)*) f(g) (g) Ao)(fc)) 

where the second step uses the dehnition of the product on 'H®"'. Applying the trans¬ 
position r to this equation will swap j -H- /c in the double tensor products. However, note 
that (^ 1 ^ fi,ij){k) = A"'®"'A"'/j = A"'Vj by [28, Lem. 1.1.9(2)]. Since "H is cocommutative, 
we have = n A^ for any A > 0 and any permutation vr: T-Ln —^ 'Hn of the tensor 

factors. In the special case of tt = r on 'H „2 = ?f®", this yields 7 -A"®”A" = A”®”A"', 
so we can keep the argument A”®"A"/j. Since r is furthermore a morphism of algebras, 
M*M* f comes to 

n V”®“ (^(gl (g) (A”®”A7i) = n (® V” (^® A”j A”/,, 

= n (® (E At)) A”/i = n ® (MM)’, km, 

where we have used (15) in the last but one step. Comparing this with (19) and using 
again the dehnition of the product on 'H®”', the claim follows. □ 

Let us give some important examples of scaled bialgebras, which will also turn out to be 
admissible coefhcient algebras for building a suitable ring of Rota-Baxter operators. 

Example 3.4. The prototypical example of a scaled bialgebra is the polynomial ring K[x\. 
Its coproduct is given by A(x) = a; 0 l-|-l 0 a;, which implies A(x") = iT) ® 
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on the canonical -ft"-basis of K[x\. In other words, we have Ap{x) = p{x + y) under the 
natural isomorphism K[x] 0 K[x] = K[x,y]. This is a well-known example of a Hopf 
algebra, sometimes also called the binomial bialgebra [17, Ex. 2.2.3.2], We use the scaling 
action x K[x\ K[x\ given by the substitutions X*f{x) = f{Xx). It remains to 
check that the action A eA A* respects the convolution product. Indeed, on the K-basis 
element x™' we have 

m , X 

(A* = V(A*®/x*) A(x™) = (’T) 

k=0 ^ ^ 

m f \ m , ^ 

= E 7 ) ® = E 7 (Ax)''(ax)"“-* = (Ax + fix)” 

k=0 ^ ' k=0 ^ ' 

= (A + /i)*x"^, 

as is required for a scaling. Under the aforementioned isomorphism 77[x]®” = K[xi ,..., x„], 
Dehnition 3.2 yields now the induced matrix action M.^[K)xK[xi ,..., x„] —)■ K[xi ,..., x„] 
dehned via (18) by M*{x\^ ■ ■ ■ x^) = for the matrix M G M.n{K) 

and exponent vector (/ci,..., kn) G N"". Since A"'/(x) = f{xi + ■ ■ ■ + x„), the row action 


^n) 

{MaxI H- MinXnfX 

induces the overall action M*{x\^ ■ ■ -x^ ) = Hi ('Ylj with the intended effect of a 

linear substitution x, i-A Yhj MijXj. In other words, the column (xi,..., x„) gets multiplied 
on the left by the matrix M G Wl„(77). 

Similar considerations apply to the larger ring of exponential polynomials K[x, and 
its variants (e.g. restricting the exponents to submonoids of K, like replacing by e^^). 
Since the coproduct is an algebra morphism, it is sufficient to dehne it on the algebra 
generators x and Clearly, K[x] should be a sub-bialgebra, so A coincides on x. For the 
exponential generators, one puts Ae“^ = so these are group-like elements unlike 

the primitive element x. Under the isomorphism iC[x, = iC[xi, • • •) 

we have again A"/(x) = /(xi,..., x„) for any / G K[x, e^^]. Moroever, it is easy to check 
that (A -I- p)* = A* 4-p*, so we have a scaled bialgebra. As in the case of the polynomials, 
one sees that M* acts via the linear substitution x* i—)■ MijXj. 

3.2. Integration in Scaled Bialgebras. Let us now turn to the interaction between the 
substitution structure and Rota-Baxter operators. Here and henceforth we shall suppress 
the unit 1 and the product V of a bialgebra (Ti, V, 1, A, I). 


M*,x’^^ = {M*,0---0 Ml) (xi -t- • • • + Xn)’^^ 

ki 

I'l} ..., f 
ki 


E 

ll~\ - \~ln — 


(M* ® . 0 M*,) {x\ 


E 

- \-ln — 


/l, . . . , 


(MiiXi)'l • • • {MinXnY" = 
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Definition 3.5. We call ("H, A, J, J) a Rota-Baxter bialgebra^ over K if [Ri, J) is an ordi¬ 
nary Rota-Baxter algebra over K and ("H, A, J) is a bialgebra such that e := 1 o J is the 
projector associated to Ini ^ R K = Ri and the horizontal substitution rule in the form 

(20) A J = ( J (g) id) A (e (g) id) A J 

is satisfied. We call {'H, A, J, J) a scaled Rota-Baxter bialgebra if it is further endowed with 
a scaling x A i—)■ A* subject to the diagonal substitution rule JA* = A“^A*J. 

The condition that e be the projector associated to Im J -(-A' = RL may also be expressed 
by 1 o J o J = o and 1 o J o 1 = 1, and this implies the direct sum. Hence we may also 
dehne a Rota-Baxter bialgebra as a a bialgebra with a Rota-Baxter operator that satisfies 
these two conditions along with (20). 

As remarked above, a scaled bialgebra has the antipode S := (—1)*, which in the case 
of a scaled Rota-Baxter algebra satisfies JS -|- Sj = 0 by the diagonal substitution rule. 
This betrays the oriented nature of this notion of integration: In typical cases, like those 
described in Example 3.4, one has S': /(x) i-A /(—x), so the integral picks up a sign under 
reflection (see also the remark at the end of Example 2.5). 

Since the bialgebra {RL, A, J) is cocommutative by hypothesis, the horizontal substitution 
rule is also valid in its symmetric variant A J = (id (g) J) A -|- (id (g) e) A J. Using (17) we 
can write both versions in the form 

(21) AJ/ = f A/ + (J/)(») = f A/ + (!f)(x) 

for functions / G "Hi C S^iRL). Here we have used the notation h{x) = h®l and hiy) = l®h 
for the two embeddings Ri RL 2 - Moreover, the horizontal substitution rule (20) can be 
iterated as follows. 

Lemma 3.6. If {RL, A, J, J) is a Rota-Baxter bialgebra over K, we have 

( 22 ) A"J = (J id®(”-^)) A" + (e ® id®^”-^)) A"J 

for any n > 0. Using the operations of Sj{RL), this may be written as A" J = -|- E^, A". 

Proof. Note that the case n = 1 is trivial since A^ = id and eJ = 0. Hence we use 
induction over n with the base case n = 2 for A^ = A, given by the hypothesis (20). Hence 
assume (22) for a fixed n > 2; we show it for n -|- 1. Using the dehnition of the base 

case and then the induction hypothesis yields 

A"+7 = (A” 0 id)(J 0 id) A + (A*" id)(E 0 id) AJ = (((J 0 id®^'^-^)) A^) 0 id) A 

-L (((e 0 id®(’^-i)) A’^J) 0 id) A + (A” 0 id)(E 0 id) AJ. 

^Note that this notion is distinct from the Rota-Baxter coalgebra in sense of [24]. There may be 
interesting relations between the two concepts but this investigation will have to wait for future work. At 
this point, let us just mention that the standard integral on polynomials is a Rota-Baxter bialgebra (see 
Example 3.8 below) but not a Rota-Baxter coalgebra in the sense of [24]. Note also that the axiom (20) 
involves the counit, unlike the axiom in [24]. 
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Using (17) and the alternative recursion for A"" gives (E(8)id'^^"'“^^) A"" = (l(8)id®*'"'“^^) A"“^ 
and (A"" (g) id)(E (g) id) Aj = (g) J. Combining this with the dehnition of A”+^ in the 

hrst summand we get 

A'^+7 = (J 0 id®'^) A"+^ + (^((1 0 id®(”-^^) A"-^) 0 id) A + 0 J. 

Using (20), its middle summand is 

(((1 0 id®(”-i)) A"-^) 0 id) (J 0 id) A = (^((1 0 id®(”-i)) A”-i) 0 id) (AJ - (e 0 id) AJ) 
= 1 0 A"J - (^((1 0 id®(”-^)) A’^-^e) 0 id) aj = 1 0 A”J - 1®*" 0 J, 

where in the last step we have used A^'^^e = and (17). Substituting this into J, 

it remains only to check that (e 0 id®"") A”+^J = 10 A”J. Applying once again (17), 
and (16) with mi = 2, m 2 = n — 1, one sees immediately that (E0id®"') A"'+^ = 10A"'. □ 

Theorem 3.7. Let ("H, A, J, J) be a scaled Rota-Baxter bialgebra. Then we obtain a Rota- 
Baxer hierarchy (^3(77), J'^")neiN with the operators = id®*^""”^^ 0 J. 

Proof. It is easy to see that the J^'" are commuting Rota-Baxter operators on SjifH). Hence 
let us now check conditions (a) to (d) of Dehnition 2.3: 

(a) Let n < m. From the dehnition of the Rota-Baxter operators, c 77 ®”^ 

holds. For showing j^"j(4m{K)* = Aim{K)* , choose any 

M = ^ with M e Mr{K) 

and / G S){T-i). We may choose k > m + r such that / G 77®*^ with / = /i 0 • • • 0 /fc- 
By (18) we have 

M*f = (et/l) • • • (eUm) (M;+i,. fm+l) • • • fm+r) 

with the rows = (0 Mi,),..., M^+r,. = (0 M^,) each having m leading zeroes. 

From (17) it follows that (elfi) ■ ■ ■ (e'^fm) = /i ®® /m G 77®*'. Unfolding the dehnition 
of Mfi+j (7 = 1) ■ • ■)''") using coassociativity (16) we have 

MUifm+i = (e®™ 0 M®) A"^+7^+, = (id®”^ 0 M®) (e®”^A™ 0 A^) Afm+j 
= e(7+,-(i))1®-0M;./^+,.(2), 

where the last step uses the relation = E{g) G 7F C 77m for g eRL, which follows 

from (17). Altogether we obtain 

= ( n E(/m+y(l)))/l 0 • • • 0 /m 0 (8) M*Jm+j,{2), 

V j=l / j=i 

from which it is clear that J^"M*/ = M*J'^'‘/ since the Rota-Baxter operator J'^" ahects 
only the prehx /i 0 • • • 0 /m- 
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(b) We must show that is an ordinary Rota-Baxter algebra over "Hn-i with 

evaluation E*. Computing the tensor kernel via [12, Prop. 2.17] we have 

Ker = (Ker id®(’^-b) + Hn-i O Ker J = 0 «) 77 + Hn-i 0 

since J is injective; we conclude that is injective as well. Since (77, J) is ordinary we 
have also 77 = IP -j- Im J, which implies 

Hn = Un-i ®n = (77„_1 (8) K) + (77„_1 ® Im J) = 77„_i + Im J"" 

so (77n, is indeed ordinary over T-Ln-i- The projector along the decomposition above is 

clearly ® e, which sends fi®---® fn-i <8 fn to e(/„) /i <8) • • • (8) fn-i G 77n-i C 77n. 

This agrees with 

E:{fl (8 • • • ® /„-l ® fn) = (et/l) • • • K_Jn-l) (07n) 

since (e^i) • • • {e*^_Jn-i) = /i<8)- • again by (17) and 07„ = = e(/„) G K 

as in Item (a) of this proof. 

(c) If r = (7 j) is any transposition, the property follows directly from the 

dehnition of the Rota-Baxter operators J*". 

(d) Finally, we have to show the three instances of the substitution rule. Among these, 
the diagonal substitution rule (8) is an immediate consequence of the diagonal substitution 
rule in (77, A, 1, J) as specihed in Dehnition 3.5. 

For the horizontal substitution rule (6) note that Tx{ei) = (ei -1- 62) © 62 so that 

Tx{ei)*{f ®g) = (ei + 62)*/ • e^g = (A/)(l ® g) = /(p © f( 2 )g, 
using again (17) for the second factor. Since ^ g) = (J/) © g, using (20) yields 
n(ei)’U(/®9) = (Al/)(1® j) = (UA/ + 1® J/)(l®9) 

= 1/(1) ® /(2)g + 1 0 (J/)g = f Tjeiru 0 9) + Gn(ei)*f (/ 0 g), 
where the last summand comes from observing that E*Tx{ei)* = ([j j)* and hence 
E*Tx{ei)*{jf) g = (e^/) (es^) = (E©id) Aj/ • (e © id) A^ = 1 © {Jf)g, 
where (17) was employed for the last step. 

It remains to show the vertical substitution rule (7), which will need a bit more effort. 
By Lemma 2.4, we may assume that v' G {0,1}"'“^. Let us write tc = Ci + n G {0,1}”“^, 
and let A C the set containing 1 and all indices i with Wi-i = 1. By the dehnition 
of and it suffices to show (7) on 77^. Hence let / = /i ©•••©/„ G 77„ be arbitrary 
but hxed. We can write this as with p = /;(©••• © f{ (j = 0,1) and 

/° = I ^ if f G A, ^ ifi if i G A, 

* I fi otherwise, * | 1 otherwise. 

Observe that /° = 1 © 1 © • • • commutes with and Using (17) we obtain also 

(23) T,r(w)*F/= T,j,(u;)7j/i) ©/2 © • • • ©/n =/°et(J/i) U (ei + d)* fi 

*ga\{i} 
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and hence L^^w)* f = A similar calculation shows that /° commutes 

also with Ly{v')*, L~^{v')* and L^^ei)*. Consequently we can pull out /° both from 
the left and right-hand side of (7), and we may thus assume without loss of generality 
that = 1 G l-Ln- Let k be the cardinality of A. If r G Sn is any permutation sending A 
to {1,, k}, conjugation by t* will reduce (7) to the case / G Tik C Tin- Hence we may 
also assume that k = n and A = {1,...,?7,}. Then (23) shows that 

= ei(J/i) n(ei + ei)*/i 

i>l 

Using again (17) we have ej(J/i) = (J/i) ® 1 ® • • • 0 1 G From the dehnition of the 
matrix action we have A = (2 i)*{ei + e*)*, which implies 

rL^iwrrf =r (d/i) n kw ® <8> kA = (j (j/i) n kw) ® (S) km 

^ i>l 2>1 i>l i>l 

= ((J/i) (J n km)) ^ (8)/*,(2) - j(/i J n km) km 

^ i>l ^ ^ 2>1 ' i>l 

(24) = (j/i)r(n A(i)® (S)A(2))-r/i((jn A(i))® <s>km) 

2>1 i>l ' ^ 2>1 i>l ^ 


where in the last but one step the Rota-Baxter axiom for "Hi was applied, while the last 
step uses the embedding "Hi C Turning now to the right-hand side of (7), we obtain 


and hence by (20) also 


LyW)*f = /i ® /2 n km ® <8 km 

i>2 i>2 


LMrrrkk'n = LMrUsk) »(j/2 n km) ® ® kA 

^ i>2 i>2 ^ 

= (J/i) (^(J/2 n A(i)) <S>km) 

= fi) r {k,m n fi,mm h,m H ft,mm ® 8 A( 2 )) 

^ i>2 i>2 i>2 ' 

+ (J/i) ® (J/2 n /l(i)) ® 8 fi,m) ■ 

i>2 i>2 ' 


Using the antipode S = (—1)*, we obtain further 


(25) - (J/i) J (/ 2 ,(i) n kmm ® km 11 fi,mm ® ® f^mm 

^ i>2 i>2 i>2 

(26) -d (J/i) <8 (J/2 n /*,(!)) n *S'/i,(2)(i) <8 8 fi,mm 

i>2 i>2 i>2 


Let US hrst study the left summand of (25), which we claim to be equal to the hrst term 
of (24). For this, it is sufficient to show that 

(27) n kmm ® n fi,mm ^kmm 8 kmm = H km ® ® 8 km- 

i>2 i>2 i>2 i>2 i>2 
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From (16) we have (A 0 A)A = A"^ = (id 0 A 0 id)A^, which may be used to rewrite the 
left-hand side as 

n(3i)*(id0 V 0 id)(A 0 S' 0 id) /*,(!) 0 /i,(2)(i) ® /i,(2)(2) 

i>2 

= n(3 *)*(id 0 V 0 id)(id®^ 0 S' 0 id)(A 0 id)(id 0 A)A/j, 

i>2 

= n (3 i)*(id 0 V(id 0 S') 0 id)(id 0 A 0 id)A^/j = H (3 f)*(id 0 e 0 id)A^/j 

i>2 i>2 

= n(3*)*(id0 E0 id)(A 0 id)A/i = 11(3*)*/*,(i) ® 1 Am 

i>2 i>2 

where the third step applies the identity V(id0S')A = e, which is true because ("H, J, A, S') 
is a Hopf algebra. The last expression follows from (17) and is clearly equal to the right- 
hand side of (27). 

We determine now the second term (26) of the commutator on the right-hand side of (7). 
To start with, we compute 

L,(e,)*f = /, a(J /2 n km) ® ® kmt 

i>2 i>2 

= /i J'^(/ 2 ,(i) n /j,(i)(i) ® h,m n fi,mm ® ^ Am') + /i ® (J /2 n fi,m) ® <8) Am^ 

^ i>2 i>2 i>2 ' i>2 i>2 

using once again (20); hence the required commutator term *JX(ei)*J%(n')*/ 

is given by 


(J/lJ/ 2 ,(l) n A(l)(l)) / 2 ,( 2 ) n A( 1 )( 2 ) 5'/i,(2)(l) (8 (S) fi,mm 

i>2 i>2 i>2 

+ (j/i) ® (J /2 n km) n sf. ,( 2 )( 1 ) «) 8 Amm- 

i>2 i>2 i>2 

We notice that the second summand cancels with the second summand of the first com¬ 
mutator term (25). Therefore it remains to prove 


(J/iJ/ 2 ,(i) n Amm) A,m H Amm ^Amm 8 Amm 

i>2 i>2 i>2 

= rA(lA,m n Amm ^ A,m H Amm ^Amm ® 8 Amm) 

^ i>2 i>2 i>2 ' 

= r/i((j n //.(!)) 8 A( 2 )), 

^ i>l i>l ^ 

where the first equality uses only the definition of J^. But this clearly follows from 


/ 2 ,(i) n Amm ® A,m H fi,mm ^Amm ® 8 ft,mm - H fkm 8 Am^ 

i>2 i>2 i>2 i>l i>l 

which is a trivial consequence of the identity (27) that we have shown above. 

This concludes the proof that J'^")nGN is a Rota-Baxter hierarchy. □ 
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Example 3.8. Let us now make sure that the scaled bialgebras of our standard model 
(Example 3.4) are in fact scaled Rota-Baxter bialgebras under the natural choice of Rota- 
Baxter operator. We do this for the exponential polynomials K[x, e^^], which includes the 
plain polynomials K[x\ as a Rota-Baxter subhierarchy. 

The algebraic dehnition of the Rota-Baxter operator J on K[x^ can be given either 
in recursive or in summation form; the latter is more amenable for our purposes. Hence 
we have for A; G N, a G JL the formulae 



OiX 



(-l)'=+ifc! 


k 



k—i ax 

Q,i+1 


i=0 


k+1 ’ 


(« 7 ^ 0 ), 


where k- = k\/{k — i)\ denotes the falling factorial. For verifying that (iL[a;, e^*], J) is a 
scaled Rota-Baxter bialgebra we have to check that it satishes the diagonal substitution 
rule JA* = (A 7 ^ 0) and the horizontal substitution rule (20). The former is 

immediate, so let us verify the latter in the form (21). For the JL-basis vector / = 
we compute 


AJ/ 

fA/ 


(-l)fc+ifc! 


k—i 




EEer) 

i=0 j=0 


k-i\ (-l)"fc- a{x+y) 


A»+i 


■ x-'y 




(~1)" rf.l-i„ik-lpa{x+y) 

„i+i -rye , 


1=0 


1=0 i=0 


k 

(_l)fc+ifc! {-ly+^kL ^,k-l„ay 

—-— y e 

1=0 


We observe that the constant terms of AJ/ cancels with the one of {jf){y), likewise the 
hrst term of J^A/ with the second term of (J/)(i/); thus it remains to show the two double 
sums equal. But this follows immediately by transforming the index i of the outer sum 
in Aj/ to I = i + j and then swapping the summations. The verihcation for the iL-basis 
vector / = is similar but simpler. 

Since every ordinary Rota-Baxter algebra {Q,j) contains K[x], its induced hierarchy 
contains the ascending algebra of polynomial rings /L[xi,... ,x„]. These are closed under 
all integrators and linear substitutions, so one can always choose iL[xi, X 2 ,... ] as the 
simplest coefficient domain. Another natural choice for the classical example C'°°(]R°°) is 
given by the exponential polynomials discussed in Example 3.8 above. Both of these are 
instances of admissible coefficient domains in the following sense. 


Definition 3.9. A Rota-Baxter hierarchy is called separated if it is of the form S){'H) for 
some scaled Rota-Baxter bialgebra Ti. If is a hxed Rota-Baxter hierarchy, an admissible 
coefficient domain for is a separated Rota-Baxter subhierarchy Q < 
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As remarked above, the minimal choice is to take Q = S){K[x]), which is an admissible 
coefficient algebra for any Rota-Baxter hierarchy {iF, We may view it as an analog 

of the prime held in any given held. 

4. The Ring of Partial Integral Operators 

We proceed now to the task of setting up a ring of partial integral operators and substi¬ 
tutions acting on a given Rota-Baxter hierarchy, and taking coefficients in any admissible 
coefficient domain. We will do this in two steps: hrst identifying hrst an ideal of suitable 
operator relations, then constructing the operator ring as a quotient algebra of the free 
algebra modulo the relation ideal. 

Note the parallel development for rings of differential operators, where the n-th Weyl 
algebra An{K) is built as a quotient of K[xi, ..., di,... ,dn\ modulo [di, xf = 5ij. In 
the case of a single Rota-Baxter operator (without linear substitutions), the corresponding 
operator ring has been studied [19] in connection with representations of Rota-Baxter 
algebras. 

4.1. Crucial Operator Relations. Logically speaking, all the operator relations that 
we shall now compile are in fact consequences of the axioms of Rota-Baxter hierarchies 
(and admissible coefficient algebras). But we need the operator relations in this special 
form since it facilitates the algorithmic treatment of partial integral operators and linear 
substitutions. Ultimately, we hope to have a noncommutative Grohner basis for the relation 
ideal (see the end of Section 4.2). 

First we need to clarify some notational conventions. Let (J^„, J'^")neN be a hxed Rota- 
Baxter hierarchy over K. For avoiding confusion with multiplication operators, we will from 
now on apply the alternative notation /[M] for the action M*{f). If M = [i j) is the matrix 
of a transposition, we write f[i j] for f[M], and parentheses are also dropped for embedded 
row vectors M. We may identify all elements f E iF with their induced multiplication 
operators iF ^ iF,u fu. Thus the operator M* f denotes the composition u i-A M*{fu) 
rather than the action f[M]. The basic operator relation M* f = f[M] M* describes how 
substitutions interact with multiplication operators. Note that /[M][M] = /[MM]. 

Let Q be an admissible coefficient domain for Then we have Q = for a 

scaled Rota-Baxter bialgebra Fi. Since the latter is determined by the coefficient domain 
as the hrst tensor grade FL = Qi, we may omit reference to Fi altogether. If / € we 
shall write /(xj) G as an intuitive shorthand for (1 i)* f. 

The crucial property of an admissible coefficient domain is that it provides tensor ex¬ 
pansion all multivariate functions g E iF as per 

r 

(28) 9 = 

/t=i 

for suitable ... ,9n,fi G IFi, where we shall consistently use Greek indices for the com¬ 
ponents. For dehniteness, we hx an ordered JF-basis (6*) for IFi, and we use the graded 
lexicographic ordering for the induced JF-basis on iFn- Then each summand above has the 
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form 0 • • • (g)so that we may choose and for j > 1. 

Consequently, each G is uniquely described by the n component sequences 


Qii] ■■= (1 if fe',1, • • •, 9j,r) e Glj) (j = 1,..., n), 

where the number r G N depends of course on the element g G J^n'-, its minimial value is 
the tensor rank of g. Note that we recover (28) by setting gj^fj_ = {g{j})^ and replacing 
the tensor product by the product of Q. For a C {1,..., n} we will also use the extended 
sequence notation 

g{a} := Y[g{j} G 

j&a 

where the product on the right-hand side is taken componentwise. If a = (ii,... ,ir) we 
denote this briefly by g{ii ,..., v}- For complements we write g{ay := g{a'}, where a' is 
to be understood as {1 ,..., n} \ a. 

We will employ component sequences in conjunction with the following variant of the 
Einstein summation convention-. Unless stated otherwise, summation is implied over all 
Greek indices on a component sequence (over the index range of the component sequence).^ 
For example, the expansion (28) can be written as g = g{l}ng{l}y Similarly, we have 
equivalent expansions like g = g{l}^g{l, 2y^g{2}^. The practical value of this convention 
will become apparent when component sequences are used in noncommutative operator 
expressions, as in the following result generalizing Lemma 2.7. Here and henceforth we 
employ the following notation for sequence substitution: If 7 G is an extended oper¬ 
ator sequence and M G M.n{K) any substitution, we write 7[M] for the sequence with 
components 7 i[M],..., 7r-[M]. In the sequel, indices i, j range over {!,..., n}. 


Lemma 4.1. Let (.F„, he a Rota-Baxter hierarchy over a field K, and Q an ad¬ 

missible coefficient domain for T. Let M G Al„(iF) have j-th column e* and g & Qi- 
Then 

M- = 7; (1 - £•) M-p 7[i jl 

for i = Hi}', 'r = 9{j} 9 = .... -Mij-i, 1, -Mij+u ..., -MiJ. 


Proof. It suffices to consider i = j = 1 since the general case follows then by multiplying 
with (1 j)* from the left and (1 i)* from the right. Hence write M = (g with v G K^-i 
and A G Ain-fiK). Note that M = (oa)(o/) with / ;= so that 


rg{xyM* = rg{i^,yao^y = 
= 7 ; (1 - El) M*f ^7^, 




^Incidentally, there is a close relation to sumless Sweedler notation: Writing 7 := G JT, for the 
(iterated) coproduct, one obtains the relation A"g = 7{l}/i ■ • • = g{i) ■ ■ ■ g{n) so that one could 

identify the “formal symbol” with the component sequence '){i}. 
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where the third equality follows from Lemma 2.8, the fourth employs tensor expansion and 
Item (c) of Lemma 2.7, the hfth follows from 7 G and Item (a) of Dehnition 2.3, the 
sixth from 7 ' G and straightness. □ 


Lemma 4.1 allows us to put arbitrary one-dimensional integrators into nor¬ 

mal form. All we have to do is to apply one sweep of Gaussian elimination to M for 
creating zeroes underneath the hrst pivot element in the j-th column of M. If the pivot 
is Mij so that My = ■ ■ • = = 0 , one uses the eliminant 


/I 0 \ 


(29) 


U{-1) = 


1 0 
1 

h+i 1 


V -In 1/ 


where I = (/j+i ,... ,ln) ^ Kn-i with := M^j/Mij. Applying this elimination yields 


(30) M:=Li{-l)-M 


with MijCi in its j-th column. 


/ Mn ■ ■ 

■ 0 ■ 

■ \ 

■ 

■ 0 • 

• Mi_i^n 

M,i • 

Mij 

■ Min 


• 0 • 

Mi^l^n 

\Mnl ■ ■ 

■ 0 ■ 

■ Mnn / 


Proposition 4.2 (Normalization of One-Dimensional Integrators). 

Let {J^n, be a Rota-Baxter hierarchy over a field K, and let Q be an admissible 

coefficient domain for T. Then for M G A4n{K) and g E Gi we have 


(31) 




M-,' - 7)M-r‘7K,i/Mj„ ii(0’ "/M OQ. 

otherwise, 


where i = min{/c | M^j 7 ^ 0}, and M G AiniK), I G K^-i as in (29), (30) if i ^ 00 . 
Moreover, we set 7 ' = , 7 = g{j} and g = g[-^, • • •, 1, ■■■, “ wi ■ 


Proof. Assume hrst i = 00 . Again it sufhces to consider j = 1 since the general case 
follows by multiplying with (1 j)* from the left. We have to show gM* f = g^M* f 
for all f E A. We may assume f E Ak for some k > n, and we may also view M G Aik{K). 
Setting r = (1 /c), we obtain a matrix Mr G A4.k{K) with /c-th column zero, and Item (e) 
of Lemma 2.7 implies that (Mr)* f = T*M*f E J^[k)- Then we obtain ^^^{T*g){T*M*f) = 
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{T*M*f) from Item (c) of the same lemma. Multiplying this relation by r* from 

the left yields the desired identity. 

Now assume the minimum exists. We use the decomposition (30) to obtain 

M* = 

such that the j-th column of di^i/MijM is e^. Applying Lemma 4.1 and ( 8 ) yields 
T'j(ij) Af = pS 

= (1 - E;M.uM„MrpiiipdiMp.{ir 
= 7 ; (1 - HP 

= A'U' 7; (1 - E]){d,_,,M.pirdlM.,Plldi,i,M.,l HP 
= M-d pi - E;)M-p^[d,,iiM,PHiy 

which is indeed (31) for z 7 ^ cxo. □ 

Before we proceed to the ordering of one-dimensional integrators, let us note the following 
natural commutativity result for integrals and evalutions along distinct axes. 

Corollary 4.3. For i ^ j the operators and E* commute. 

Proof. Applying Lemma 4.1 with M = Ej and g = 1 yields ^^"E* = {1 — E*) E* which 
implies commutativity because of EiEj = EjEi and (f) of Lemma 2.7. □ 

By Proposition 4.2, the basic building blocks are one-dimensional inegrators of the 
form ^^'‘h[xi)Li{v)* along a coordinate axis Xi. We shall now investigate how two such 
integrators interact with each other, starting with the case of different coordinate axes. 
It turns out that we can always reorder them. In this paper, we choose to put them in 
ascending order. 


Proposition 4.4 (Ordering of One-Dimensional Integrators). 

Let {En, J'^")neN be a Rota-Baxter hierarchy over a field K, and Q be an admissible coeffi¬ 
cient domain for E. Then for i < j, vectors v G w G and functions g,h ^ Qi 

we have 


(32) ph{xj)Lj{wYpg{xi)Li{vy = (1 - Eypi^LYv')*pri^Lj{w)Y 

where we set v' = Ljf-{w) v = (uj+i, • • • , Vj, Vj^i — Wj^iVj ,..., — WnVjY G K'^~\ as well 

as Tj' = g{xi) h{i}, rj = h{j} and h = h[ej — VjCi] G G(ij). 


Proof. First note that Lemma 4.1 with g = 1 yields j^Ej{w)* = Lj{w)*p — (Lj{w) Epp, 
but the last summand vanishes by (f) of Lemma 2.7 since Lj{w) Ei = Ei Lj{w). Hence 

h{xj)LYw)* p g{xi)Li{vY = P'h{xj) p lywY g{xi)Li{vY 

= Pa(HrVHi)EP)‘Lj{w)’ = (1 - Ej) HP'pd'O}^ UW* 

= (1 - E;)pgP) hPMvj-pMALAP’ 
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where the second step follows from (4) and Item (c) of Lemma 2.7, from Lj[w)_,i = li and 
the commutativity of Rota-Baxter operators, the third from Lemma 4.1 with h = Li{v')*^h 
and Li{v')j, = —+ = Cj — vjej, the last step from Corollary 4.3 and straightness. □ 


Lemma 4.5. Let (Jm, J'^")neiN be a Rota-Baxter hierarehy over a field K, and Q an ad¬ 
missible coefficient domain for T. Then for i ^ N, h E Qi and w G \ {0} we have 


(33) r'Hxi) L,(wyf“ 




P‘U(wy)Pr,,Lt{w'r 


as an identity on Tn, where k = min{j > i \wj t)} and rj' = h{n + 1 }, = h{k} with 

h ■ /tc/j]/tc/j G 

Moreover, for w = (tCj+i,..., Wnfi we have defined 

W = WkCk-i G K^~\ w' = (Wfc+i, . . . , Wn)/Wk G 


as auxiliary vectors. 


Proof. It suffices to assume f = 1; the general case follows by conjugation with r* = [i 1)* 
since rLkiw'ffi = Lk{w') and tLi{w)t = Lfiw). We write J := /„©efe as an abbreviation. 
Using the relation L^iw)* = wlL^^Ck-i + [1/wk)*, we obtain 


rh{x)L,,{wYj'" = w^^wl rHdi,i/wY LY^k-i + (0, 

= wi Lpw'yr (i.(e»-.)T - f ‘A lywy (^y 

= Lp(«.')V*(«>U4et-i)T - u,lS‘L4et.y)pkU(w'y (Jj) 

= Ly\w'y.r ffiwy f - rLywy)phwiLyw'y 


* 


where the hrst equality uses ( 8 ) twice, the second employs (13), the third the commu¬ 
tation of tCfceii + 622 + • • - Grin with {w') aud J, the fourth applies again ( 8 ) together 
with L 3 ;(eA:-i)* = Lx{w)*wl, and the last wlLk{w')*{l/wk)* = Lk{w')*. Finally, (33) fol¬ 
lows by substituting the expansion of h, noting that the components rj' = h{n+l} G 
commute with and L^iw)* and that J* acts on them as {k n + 1)*. □ 


For stating the last operator relation, we need the following technical result for composing 
two substitutions. A general result could be formulated but for our present purposes the 
following special case will be sufficient, where the hrst substitution is Xy i-A axk + bxi, and 
the second is Xk eA cXi - 1 - dxj. 

Lemma 4.6. Let {Qn, J'^")neN be a separated Rota-Baxter hierarchy over K. Fix h G Qi, 
indices i < j < k < I, and scalars a,b,c,d G K. Assume h[aek + hei] = u is expanded 
into u = ri'^ryj_ with p = u{k} G and p' = u{l} G Setting p^ = u^{xk), expand 


MULTIVARIABLE INTEGRATION AND LINEAR SUBSTITUTIONS 


35 


also u^[cei + dcj] = with 7 ^ = e and 7 ^ = u{i} e %nto 

where summation is only over n. Then in we have the standard tensor expasion 

(34) h[acei + adej + bei] = 7 ^ 7 ^,^ 7 ^^, 

with the sum over (//, G {1,..., Si} x • • • x {1,..., = {1,..., Si + • • • + s^}. 

Proof. Choosing n > I, write M = ack+bei G A^„(iC) and M = Ik-i(B{cei+dej) G M.n{.K) 
in the sense of the embedding. Operating M* on h[M] = rj'^rjfj, yields 

h[acei + adcj + bei] = h[MM] = 

since M* is an algebra homomorphism which leaves rj'^ G Q(i) invariant (this follows easily 
from straightness and M;, = ei). But = r]^\ik\ implies = rj^fM] since the matrices M 
and (1 k) ((ccj + dcj) © In-i) have the same k-th. row (again the conclusion follows by 
straightness). Now (34) follows by substituting the tensor expansion of u^. □ 

We can now state the last important relation, which will enable us to put one-dimensional 
integrators into strictly ascending order since two adjacent integrators along the same 
coordinate axes can be reduces as the following propostion makes precise. 


Proposition 4.7 (Coalescence of One-Dimensional Integrators). 

Let (iFn, J'^")nGN be a Rota-Baxter hierarchy over a field K, and Q an admissible coefficient 
domain for F. Then for i G N, and v,w E iC"”* with w ^ 0 we have on Fn the identity 

^'"'h{xi) Li{w)*pg{xi) Lfiv)* = Ll^{w')*:x"[k n + l]^ x 

^ ^ ^ pg{xi)x'^,Li{v')* - pg{xi)x'^,Li{v'+ w)*'^J'"^Xf^Lk{w'y, 

where we have set 

h = h[-{v'^./Wk)ei + Cj/Wk - Cn+l/Wk] 

v' = V = {Vi+i, ...,Vk, Vk+i - w'k+^Vk, ...,Vn- w'^VkY G 

with component vectors x" = h{n + 1}, x' = h{i}, x = h{k}, and w', w from Lemma f.B. 


Proof. Multiplying (33) by g{xi) LfivY from the right, the right-hand side of (35) becomes 


Lk\w'Yi[k n+lY i^LfiwY P - pLi{wYjPVf.Lk{w'Yg{xi)Li{vY- 
Since Lk{w'Y commutes with g{xi) and (4) swaps eliminants, the right-hand factor is 

pr]f, Lk{w'Yg{xi) LfivY = g{xi) Lfiv'YLk{w'Y 

= 9{xi) ipi - El) Li{v'Ypx^,v Lk{w'Y, 


where the last expression (summation only over u) follows from Lemma 4.1 by writing 
Vii ~ and defining = ufick with component vectors = g{i}, Xr. = 9{k}- 

By Lemma 34 we can determine s G N and make the identihcations rj' = x" ^ 
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7 ' = ^ ^ 7)5 and X = 7 £ Gfk)- Comparing the above equations with (35) and noting 

that g{xi) x'fj, G G{i) commutes with Li{w)*, it is clear that it suffices now to prove that 


[Li{w)*p - pLi{wy^g{xi)x'^,El = 



- PHP*El^g{xi)x'f, 


vanishes, where we have used that g{xi)x'Pk — g{xi) x'fj, by g{xi)x'fj, ^ G{i)- But the 
parenthesis on the right-hand side above is zero by Li{w)*El = and Corollary 4.3. □ 


We have now collected all essential operator relations for building up the partial integral 
operators as a quotient ring modulo these relations. 


4.2. Construction via Quotient. The raw material for building the operator ring over 
an admissible coefficient domain Q for some Rota-Baxter hierarchy E comes from the 
three base algebras: the coefficient algebra Q] the substitutions from K[A4] in the form 
of a monoid algebra, and the commuting Rota-Baxter operators (a polynomial ring in 
inhnitley many variables). These three algebras are put together by a free product, then 
we take a suitable quotient modulo the ideal of operator relations. 

We denote the free product [10, §111.6] [16, §29] [23, p. 85] of two K-algebras A and B by 
AUkB. Rougly speaking, this consists of R'-linear combinations of alternating products 
of elements of A and R; for a short introduction see also [15, Example 3.28] where K 
is the complex held (though the construction works for general helds) or the historical 
reference [22]. One has natural isomorphisms A 11k {B BaC) = {A ILkB) 11k C so that 
we may employ the notation A 11k B 11k C. Incidentally, note that the free product is 
commutative in the sense that A 11k B = B 11k A, despite having a noncommutative 
multiplication. 

Recall that A1 := A1{K) denotes the substitution monoid over K, which contains the 
semigroup = M(K) \ {/}. Since its action on E and Q < E rs contravariant, we 
form the opposite monoid algebra K[AA]* and the associated opposite semigroup algebra 
K[Aiy*. Its basis elements will be written as M* with M G K[A1]. For avoiding confusion, 
we shall write the integrals in the prospective operator ring as An (n G N); they form the 
commutative polynomial ring K[A\ = K[An | n > 0]. After setting up an appropriate 
action 0 on the given Rota-Baxter hierarchy (R„, J^")nG]N 5 we shall have An 0 / = Pf 
and M* Q f = f[M] and f Q g = fg. 

Definition 4.8. Let {Q, J) be a separated Rota-Baxter hierarchy over a held K. Then the 
ring of partial integral operators over Q is dehned as the quotient JL-algebra 

(36) g[y := g 11kK[MY 11kK[A] / ^ 

with the congruence = generated by Table 1 for the implicit ranges i, j G N such that i < j 
and M G K[fA\, and g,h ^ g\. For all other notations and conditions, we refer to the 
Propositions 4.2, 4.4, and 4.7. 

Note that Table 1 is meant to contain various implied congruences as special instances 
of the ones given, obtained by substituting g{xk) = 1 in the hrst rule; g{xj) = 1 in the 
rules 4, 5, 6; and h{xj) = 1 or g{xi) = 1 or both in the last two rules. Moreover, n = 0 and 
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M*g ^ g[M] M* 

M*Ai ^0 if Mi. = 0 

Aj 9 {xi) = 9 {xi)Aj 

Agixj) = g{xj)Ai 

A, g(x,)M- - My G '(1 - i.( 0 * 

if i := min{/c Mkj 7^ 0} 7^ 00 

Ajg{xj)M* ^ {pg{xj))M* 

if i ;= min{/c Mkj 7^ 0} = 00 

Ajh{xj)Lj{w)*Ai g{xi)Li{v)* = (1 - E*)Air]'^Li{v')*Aj 
Aih{xi) Li{w)* Ai g{xi) Li{v)* = n+l]^ x 

VuLj{w)* 

X {L,{wr pg{xi) x’.Hv'r - pgixi) x',Hv' + («;')* 

Ajg{xj)Aj = {pg{xj))Aj - Aj{pg{xj)) 


Table 1. Relations for Partial Integral Operators and Substitutions 


hence Liiy) = / is allowed in the last rule. These implicit rewrite rules are also important 
for the implementation. 

The definition of the free product H = Hi ]1 k ■ ■ ■ Hk yields the coproduct in the 
category of iP-algebras, with insertions Hk ^ H. This requires for each factor a unitary 
splitting H^ = K (B Hi as JP-vector spaces such that the alternating products are built 
from Hi rather than Hi, with being the empty product. In our case, the implicit rules 
will remove the units Ig & Q and I* G K[M]* and 1x[a] ^ K[A\ in favor of the new unit 
in iP C ^[J]. We use the unitary splittings 

Q = K ® U„>o Qn^ 

K[M]* = K®K[M^]*, 

K[A] =K®{pe K[A] I deg(p) > 0} 

for the three factors in the free product (36). 

A detailed description of the computational realization of Q[j] will be the subject of 
another paper. For our present purposes it suffices to observe that we can start from the free 
R-algebra generated by a basis (&i)ig]N of Q such that bo = 1, the substitutions M* G K[M] 
and the integrators {An)n>o- Then one has the following extra rules in addition to the 
explicit and implicit rules of Table 1: 

• Every occurrence of ho, I*, and 1k[A] is removed, while the empty word is treated 
as the true unit 1 G G[j]- 

• After operating on a basis function bi E Q hj multiplication, integration or substitu¬ 
tion, the resulting element of Q is again expanded with respect to the basis 

• Within monomials of A"[A], one enforces the ordering Ai > A 2 > ■ ■ ■, thus achieving 

normal forms A“ = A'^^A^^--- for «: N. In view of Table 1, reduction 

ensures G {0,1}. 
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For the subtle interplay between canonical forms and basis expansion, we refer to [37, §4.4], 
where this is discussed for ordinary integro-differential operators (but the same applies here 
mutatis mutandis). 

Proposition 4.9. Let (JTn J^")neN be a Rota-Baxter hierarchy over a field K, and let Q 
be an admissible coefficient domain for T. Then the action 0: ^[J] x —)■ induced 
by g Q f = gf, M* 0 / = f[M] and AiQ f = is well-defined. 

Proof. As indicated, one may interpret ^[J] as a quotient of the free iF-algebra generated 
by a iF-basis (fei)je]N of Q, the substitutions M* E K[M.] and the integrators {An)n>o- 
Since the implicit rules and the three extra ones listed above are obvious consequences, it 
suffices to verify Table 1 for basis functions g,gj,gk. Rule 1 is valid because K[M] acts 
on ^ C via iF-algebra homomorphisms, by hypothesis (Dehnition 2.1). Rule 2 follows 
from M = EiM and Item (f) of Lemma 2.7, Rules 3, 4 by Item (c). Rules 5, 6 correspond 
to Proposition 4.2, Rule 7 to Proposition 4.4, Rule 8 to Proposition 4.7. Finally, Rule 9 
follws since each is a Rota-Baxter algebra in its own right. □ 

We orient the congruences of Table 1 from left to right. Moreover, we agree that 
on the right-hand side all matrix products, substitutions and coefficient multiplications 
are immediately carried out. For example, the result of Rule 3 is a iF-linear combi¬ 
nation of monomials having the form Q ■ K[Ai]* ■ Ai ■ ■ K[A4\*. For showing ter¬ 

mination of the induced rewrite system (including the implicit and extra rules), we use 
the following term order on the underlying free algebra. Every monomial is segmented 
into W = tCoAjjtCiAjj • • • Ai^Wn such that the words wq, tci,..., tCn do not contain any A. 
Note that every word Wk with /c > 0 has a corresponding integration index ik, and by 
convention we set io := 0. For any subword Wk, let a{wk) denote Wk viewed as an 
element in the following Noetherian partial order on the word monoid over the alpha¬ 
bet Q l±) K[M]*. We impose no order amongst the g E Q or amongst the M* E K[M]* 
but we stipulate that g < M* whenever M ^ I] this partial order on the alphabet is 
extended to a graded lexicographic partial order on the corresponding words. Further¬ 
more, let N{wk) be the total number of nonzero matrix entries with r < ik and 

nonunit matrix entries for all M* occurring in the word Wk. If Ui is another word, 

we put now Wk < ui iff {N{w),—ik,o'{wk)) < {N{u),—ji,a{ui)) in the pure lexicographic 
sense. Finally, let U = uoAjj^uiAj^ ■ ■ ■ Aj^Um be a second monomial in the free algebra. 
Then we dehne hF < t/ iff {wn, • • •, wi, wq) < (um, ■ ■ ■ ,ui, uo) in the sense of the graded 
lexicographic order induced by the previous order given on the components. 

Theorem 4.10. Let (Q, J) be a separated Rota-Baxter hierarchy over a field K. Orienting 
the rules of Table 1 from left to right, one obtains a Noetherian reduction system. 

Proof. We note hrst that all the rules of Table 1 respect the above term order in the sense 
that all monomials on the right-hand side are smaller than the monomial on the left-hand 
side. For Rules 1-4 this is evident. Rule 5 is only applicable if M does not have the 
form Lj{v) so that N{w) > 0 for the word on the left-hand side while clearly N{u) = 0 for 
all rightmost words of the right-hands side monomials. In Rule 6 we have a drop in the 


MULTIVARIABLE INTEGRATION AND LINEAR SUBSTITUTIONS 


39 


number of integrators, so the left-hand side is greater due to our choice of graded order 
on subword sequences. Looking at Rule 7 we see that the last subword w has N{w) = 0 
before and after reduction but its negated integration index drops from —i to —j. The 
same is true for Rule 8, where the index drops from —i to —k. Finally, Rule 9 decreases 
the number of integrators and hence drops in the graded order on subword sequences. 

Next let us convince ourselves that > is a term order. If 1 denotes the empty word, we 
have 1 < W for all non-empty words. This is because the empty word has N{1) = 0 and 
integration index 0, and cr(l) is the mimimum among all a-words. Now we assume W < U 
and we must show WV < UV and VW < VU for all words V in the free algebra. But this is 
clear since multiplication by V on either side equally affects N{w) and the integration index 
of w of each segmental word w within V, and the order on the a{w) respects concatenation. 

It remains to prove that > is Noetherian. By way of contradiction, assume an inhnite 
descending chain exists. Since > is graded on subword sequences, there must then be 
an inhnite descending chain of words woAi^wiAi^ ■ ■ ■ Ai^Wn with hxed sequence length n. 
Since the N{wk) and the integration indices cannot drop inhnitely often, we obtain an 
inhnite descending chain of cr{wk) instances. But this is impossible because the order on 
the u{wk) is Noetherian. □ 

In a Noetherian reduction system, every element has at least one normal form. Each 
normal form may be taken to be a representative of the corresponding congruence class 
of the quotient algebra (36). In our case, we can go one step further: We can specify 
monomial normal forms, such that every element of the quotient algebra can be written as 
a iF-linear combination of these (congruence classes of) monomomial normal forms. This 
can be done as follows. 

Let us call a word (of the underlying free algebra) a line integrator of index i if is has 
the form Aih{xi)Li{v)* with a iF-basis element b G Qi, where again b = 1 or Li{v) = I 
means absence. A volume integrator is word of the form K = bM*Ji ■ ■ ■ Jr, where b E G 
is a basis element and M* G K[A4]* with ^ 0 for r > 0, and where Ji,..., Jr are line 
integrators of indices R < • • • < v- 

Theorem 4.11. Let {G,j) be a separated Rota-Baxter hierarchy over a field K. Then 
every element of G[J] is a K-linear combination of terms b M*Ji ■ ■ ■ Jr, where b E G is 
a basis element and M* G K[M.]* with ^ 0 for r > 0, and where Ji,...,Jr are 
integrators of indices R < ■ ■ ■ < ir- 

Proof. We have to show that every word can be reduced to a linear combination of the 
form stated in the Theorem. Using Rule 1 will ensure all subwords Wk in the canonical 
segmentation w^Ai^wiAi.^ ■ ■ ■ Ai^Wn have the form b^ Mf for some b^ E G and G K[JA\. 
By Rules 3 and 4 we can achieve that b^ G G(it,) and by Rules 5/6 that = Li^{v)* 

for some n; again n = 0 means absence. Hence all monomials can be reduced to the 
form bM*Ji ■ ■ ■ Jr, where 6 G ^ is any JL-basis elements and M* E K[JA\* is arbitrary. 
If the integration index drops between any neighboring integrators Ji, Ji+i, application of 
Rule 7 hxes this. Hence we can ensure that the sequence of integration indices is increasing. 
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For achieving a stricly increasing sequence, it remains to apply Rule 8 if tc 7 ^ 0 and Rule 9 
otherwise. Finally, Rule 2 is used to guarantee 7 ^ 0 if r > 0. □ 

In fact, it is our conjecture that the reduction system induced by Table 1 (or perhaps 
a slight variation of it) is not only Noetherian but also confluent (modulo the JF-algebra 
axioms, as usual). In other words, we have a noncommutative Grobner basis for the 
relation ideal dehned by the noncommutative polynomials corresponding to these rules. 
As a consequence, the volume integrators would be a JF-basis of G[j]. 

Conjecture 4.12. Let be a separated Rota-Baxter hierarchy over a field K. Ori¬ 

enting the rules of Table 1 from left to right, one obtains a convergent (i.e. Noetherian and 
eonfluent) reduction system. 

The confluence proof (of which we have completed signihcant parts) appears to be rather 
lengthy and laborious, at least if one follows the well-known strategy of the Diamond 
Lemma [7]. Moreover, such a proof will have a much more computational flavour than the 
current article. Therefore we plan to present the confluence result in a separate publication. 

5. Conclusion 

As pointed out earlier, our construction of G[j] is originally motivated from eomputa- 
tional algebra: We want to build up algorithmic support for calculations with integral 
operators, substitutions and function expansions. In particular, we aim to represent, com¬ 
pute and manipulate (e.g. factor) Green’s operators for suitable classes of LPDE boundary 
problems, as explained in the Introduction. This necessitates the following further tasks: 

• Adjoining derivations to form partial integro-differential operators with linear sub¬ 
stitutions. 

• For both operator rings, prove confluence to ensure unique normal forms and hence 
deciding equality (relative to deciding equality in the given coefficient domain Q). 

• Develop methods for computing and factoring LPDE Green’s operators, as for 
example in [34]. 

• Study how the classical Green’s function can be extracted from the Green’s operator 
of certain LPDE boundary problems. 

• Provide an implementation in a computer algebra package. 

We plan to address these and related items in future work. 

Apart from these rather obvious further developments, the material presented in this 
paper also contains some more intrinsic topics that might deserve a detailed investigation. 
Our notion of Rota-Baxter hierarchy {D, is a hrst attempt to provide an algebraic 

framework to study integral operators and linear substitutions. But one could also try 
laying the foundations deeper, for example by using operads with a scaling action. The 
operad of Rota-Baxter algebras [2, 3] has played a key role in understanding splitting of 
operads initiated by Today [26]. In our case, note that the operation f{xi) i-A /(xi -|- X 2 ) 
in the hierarchy D is similar to a coproduct except that it is of type T\ D 2 rather 
than ^ Di ® ^ 2 - Working on multivariate functions f{x\,... ,xfl) and using similar 
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substitutions 

(37) Aj := Ai^n- Xj H- Xj {j < i), x* H- + Xi+i, xj hA Xj+i {j > i) 

one gets a family of operations Ai^n '■ t J^n+i generalizing the coproduct of a bialgebra. 
One may check that they satisfy the commutation rule Aj Aj = Aj+i Aj for i > j. If one 
writes a\ for the addition map, one has Aj^„(/) = / o„ 2 ,i « for / G Using 

the language of PROPs [27, 29, 43], the commutation rule follows from the axioms and the 
associativity of a. In a similar vein, we can express the scaling action by = fon,i,i^* 

where X*: J^i —?■ J^i is the substitution x ha Xx. From the Ai^n and the di^x one can 
then build a Rota-Baxter hierarchy by decomposing matrices into permutation matrices, 
diagonal matrices and evaluations. 

The advantage of this approach is that it would lead naturally to a generalization 
where / j g is allowed for any g G not only for = a and g = X*. For the 
Rota-Baxter operators we would then require the general substitution rule of multivariable 
calculus (and for the derivations the general chain rule), as it is satished for C'°°(]R°°). 

Another possible generalization is more geometrically motivated: As mentioned earlier 
(Remark 2.6), for U°°(]R°°) one can interpret our formulation of the chain rule as integration 
over certain parallelepipeds. However, this is not sufficient to provide a genuine theory of 
integration over simplices. In fact, it cannot: The dehnition of simplices requires some 
notion of orientation (e.g. coming from the order structure in the case iP = R) while we 
have only assumed algebras over a general held K. Nevertheless, it seems what we are 
missing is “not much” in a certain sense. Perhaps oriented sets [6, §2.6, Expl. 9b] could be 
used for endowing the vertices of a simplex with an orientation. On a more fundamental 
level, abstract simplicial complexes might provide just the right structure. In either case, 
the simplex would have one “variable vertex” (running along a hxed 1-cell), to ensure an 
indehnite integral and hence some kind of Rota-Baxter operator. 

This in turn could be compared to the integration of differential forms over simplices, 
a species of definite integral that is crucial in de Rham cohomology (where /c-chains are 
built from fc-simplices). On the manifold R, the connection between the Rota-Baxter 
operator of indehnite integration and the dehnite integration of diherential forms is given 
by the Fundamental Theorem of Calculus. Roughly speaking, its generalization to higher 
dimensions is Stokes’ Theorem. It would be interesting to explore this connection in the 
light of Rota-Baxter hierarchies. 
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